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EXECUTIVE  SUMMARY 


Coherent,  multiple-input  multiple-output  (MIMO)  extensions  to  radar  systems  offer  a  number  of 
performance  advantages  over  more  conventional  approaches.  For  ground  moving  target  indication 
(GMTI),  the  benefits  can  include  lower  minimum  detectable  velocity  (MDV)  and  more  accurate  target 
direction  estimation.  In  part,  these  benefits  are  afforded  by  the  superior  clutter  cancellation  offered  by  the 
long,  filled  virtual  array  that  MIMO  can  yield,  while  using  relatively  few  physical  antenna  elements.  In 
addition,  the  use  of  uncorrelated  waveforms  for  coherent  MIMO  GMTI  diffuses  the  transmitted  energy 
uniformly  over  a  broad  azimuthal  coverage  region,  yielding  a  long  effective  coherent  integration  time. 
This  enhances  the  above  benefits  and  also  provides  finer  Doppler  resolution.  For  applications  in  which  the 
geometry  and  MIMO  radar  system  design  result  in  ground  clutter  occupying  a  suitably  limited  extent  of 
range  and  Doppler,  waveform  sets  such  as  time-division  multiple  access  (TDMA)  or  Doppier-division 
multiple  access  (DDMA)  can  be  used  to  realize  the  above  benefits.  However,  it  has  been  noted  that  the 
clutter  cancellation  is  less  effective  when  the  clutter  has  a  greater  range-Dopplcr  extent.  Various 
alternative  low-correlation  waveform  sets  have  failed  to  mitigate  this  problem. 

In  this  report  it  is  shown  analytically  that  the  extra  difficulty  of  clutter  cancellation  referred  to 
above  is  attributable  to  an  unavoidable  increase  in  the  number  of  MIMO  degrees  of  freedom  required  to 
cancel  the  ground  clutter,  when  the  clutter  occupies  an  extended  domain  of  range  Doppler.  Specifically, 
the  following  is  a  section-by-section  synopsis  of  the  discussion  and  results  herein. 

1.  In  the  Introduction,  time-domain  computer  simulation  results  are  presented,  showing  how  the 
signal-to-interference-plus-noise  ratio  (SINR)  for  MIMO  GMTI  is  less  than  ideal  for 
extended  range  domains  of  ground  clutter.  Specifically,  practical  waveforms  chosen  to  have 
low  cross-correlation  functions  over  extended  range  domains  are  seen  to  result  in 
significantly  worse  SINR  than  do  hypothetical  ideal  waveforms.  The  ideal  waveforms  are 
simulated  as  having  matching  autocorrelation  functions  and  zero  cross-correlation  functions. 
To  explain  this  effect,  a  pictorial  discussion  is  given,  showing  how  matched  filters  (MF)  for 
different  waveforms  effectively  weight  the  complex  clutter  reflections  from  the  ground  with 
different  correlation  functions  versus  range.  The  different  correlation  functions  result  in 
decorrelation  of  the  outputs  of  the  filters  matched  to  the  different  waveforms.  In  turn,  this 
causes  degradation  of  the  ability  to  cancel  the  clutter  returns  when  the  matched- filter  outputs 
arc  adaptively  weighted  for  that  purpose.  The  points  made  in  this  introductory  section  are 
quantified  analytically  in  succeeding  sections. 

2.  The  assumptions,  limitations,  and  applicability  of  the  ensuing  analytical  results  are  discussed. 
The  analysis  considers  cases  in  which  the  MIMO  virtual  array  is  critically  spaced  (i.e..  a/2) 
and  filled.  In  particular,  the  results  consider  only  the  clutter  on  the  “ray”  isoDoppler  to  the 


hypothesized  target.1  The  most  familiar  interpretation,  and  that  used  in  this  Executive 
Summary,  would  be  that  the  analytical  results  apply  to  a  pulsed-Doppler  paradigm,  in  which 
each  waveform  consists  of  a  sequence  of  identical  pulses  over  the  coherent  processing 
interval,  and  Doppler  filtering  (to  focus  on  the  hypothesized  target  Doppler)  is  performed 
after  matched  filtering  and  prior  to  adaptive  weighting  for  clutter  cancellation.  However,  a 
broader  interpretation  is  that  the  waveforms  are  free  to  be  of  any  form  (pulsed  or  not),  and  the 
analysis  shows  the  statistical  structure  of  the  clutter  arising  from  just  the  most  critical  region 
on  the  ground  (i.e.,  isoDoppler  to  the  target). 

3.  A  mathematical  model  is  developed  for  the  response  of  the  MIMO  receive  antennas  and  MFs 
to  the  ray  of  clutter  on  the  ground  isoDoppler  to  the  hypothesized  target.  Each  of  the  R 
receivers  is  followed  by  T  matched  filters,  corresponding  to  the  T  transmit  waveforms.  The 
RT  outputs,  sampled  at  the  range  of  the  hypothesized  target  and  Doppler  filtered  to  focus  on 
the  hypothesized  target  velocity,  are  adaptively  weighted  to  cancel  the  clutter  on  the  isoDop, 
while  maintaining  gain  on  the  hypothesized  target  azimuth  angle.  The  study  considers  the 
rank  and  eigenstructure  of  the  RT  x  RT  covariance  matrix  between  the  various  received, 
processed  outputs,  as  resulting  just  from  the  clutter  on  the  isoDop  ray  on  the  ground.  When 
more  eigenvalues  are  positive,  more  degrees  of  freedom  potentially  need  to  be  expended  to 
cancel  the  ground  clutter.  In  this  section,  it  is  shown  that  for  uncorrelated  signals,  the 
frequency  spectra  cannot  overlap  and  the  autocorrelation  functions  are  orthogonal  to  one 
another,  which  causes  decorrelation  of  the  ground  clutter  returns  at  the  outputs  of  the  MFs  for 
the  different  signals.  It  follows  that  “ideal”  waveform  sets  do  not  exist  over  an  extended 
domain  of  range.  This  last  fact  was  already  known,  but  the  extent  of  its  impact  on  the  ability 
to  cancel  clutter  was  not  fully  realized.  The  MIMO  response  to  a  point  target  is  shown  to 
have  the  familiar  Vandermonde  virtual-array  characteristic  of  MIMO,  while  the  response  to 
clutter  can  have  multiple  modes.  Also,  it  is  shown  how  single-input  multiple-output  (SIMO) 
GMTI  yields  a  clutter  covariance  having  only  one  nonzero  eigenvalue. 

4.  Under  simple  assumptions,  it  is  shown  that  the  MIMO  clutter  covariance  has  T  equal, 
positive  eigenvalues  whenever  the  T  waveforms  are  uncorrelated.  This  occurs  independent  of 
the  direction  of  the  ray  of  clutter  isoDoppler  to  the  target.  This  phenomenon  is  a  worst  case 
for  clutter  cancellation.  Also,  it  is  shown  that  in  this  case  the  receive  array  alone  performs  all 
clutter  nulling  for  MIMO,  for  each  transmitter  separately.  This  makes  MIMO  and  SIMO 
equivalent  for  clutter  cancellation,  in  this  respect. 


In  some  system  designs  with  low-velocity  platforms,  there  can  occur  hypothesized  target  velocities  for  which  no 
isoDoppler  clutter  ray  exists  on  the  ground.  Clutter  cancellation  may  not  be  required  in  such  cases.  Also,  if  there  are 
multiple  rays  isoDoppler  to  the  target  (i.e.,  Doppler  ambiguities)  within  the  azimuthal  coverage  region,  clutter 
cancellation  requires  more  degrees  of  freedom  in  order  to  place  antenna  array  nulls  in  multiple  directions.  The 
analysis  herein  considers  only  cases  in  which  a  single  isoDoppler  clutter  ray  occurs. 
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5.  For  the  case  in  which  T  -  2,  a  penalty  function  is  derived,  which  evidences  quadratic  growth 
of  the  product  of  the  highest  two  eigenvalues  of  the  MIMO  clutter  covariance  whenever  the 
two  signals  are  not  proportional.  This  result  relates  directly  to  the  size  of  the  second 
eigenvalue  of  the  covariance,  which  requires  an  extra  degree  of  freedom  for  cancellation.  The 
penalty  emphasizes  waveform  differences  at  frequencies  where  the  energy  density  spectrum 
of  either  waveform  is  high,  suggesting  the  use  of  flat-spectrum  waveforms  The  following 
statement  is  also  proven:  For  T  =  2,  the  MIMO  clutter  covariance  is  of  rank  I  or  0  for  all  <|>  if 
and  only  if  (iff)  the  two  waveforms  are  proportional.  Here,  (j)  is  defined  as  the  phase 
difference  between  adjacent  transmitters,  as  would  be  seen  on  the  ray  of  clutter  isoDoppler  to 
the  target. 

6.  For  arbitrary  T,  it  is  shown  that  the  MIMO  clutter  covariance  is  of  rank  1  or  0  for  all  <|>  iff  the 
T  x  T  waveform  cross-correlation  matrix  as  a  function  of  delay  satisfies  the  equality 
C(  r)  =  C(0)/(r),  for  all  relevant  values  oft,  for  some  scalar  (autocorrelation)  function, 
/( x ).  This  property  is  called  delay  similarity  of  the  waveform  cross-correlation  matrix. 
Using  the  above  equivalence  proves  that  the  statement  in  (5)  above  is  true  for  any  value  of  T. 
Further,  the  properties  of  analytic  functions  then  are  used  to  extend  the  statement  finally  to 
the  following  pair  of  statements:  ( 1 )  For  any  T,  if  all  the  waveforms  are  proportional  then  the 
MIMO  clutter  covariance  is  of  rank  1  or  0  for  all  <{>.  and  (2)  for  any  T,  if  any  pair  of  the  T 
waveforms  are  not  proportional,  then  the  MIMO  clutter  covariance  can  have  rank  I  or  0  only 
for  isolated  values  of  <|>. 

7.  Three  different  potential  techniques  are  examined  briefly  for  improving  the  cancellation  of 
clutter  over  extended  domains.  The  three  techniques  are  adaptively  weighted  slow-time  taps 
(i.e.,  Space-Time  Adaptive  Processing,  or  STAP),  adaptively  weighted  last-time  taps,  and 
nonconventional  reception  (i.e.,  reception  not  constrained  to  use  matched  filters  for  the  T 
waveforms).  For  uncorrelated  MIMO  waveforms,  it  is  shown  that  none  of  these  techniques 
addresses  the  basic  decorrelation  problem  of  the  MIMO  clutter  covariance.  Thus,  while  these 
techniques  could  improve  performance  somewhat,  they  are  not  expected  to  restore  near-ideal 
performance. 

8.  For  clutter  occupying  only  a  suitably  limited  domain  of  range  Doppler,  it  is  already  known 
that  near-ideal  clutter  cancellation  can  be  achieved  through  appropriate  waveform  design,  as 
noted  above.  Using  the  example  of  TDMA  waveforms,  one  can  show  here  that  it  is  possible 
to  achieve  the  delay-similarity  property  defined  in  (6)  above  when  the  clutter  occupies  only  a 
suitably  limited  domain,  which  implies  that  the  MIMO  clutter  covariance  will  have  rank  1  or 
0  for  all  <]>.  Thus,  ideal  sets  of  waveforms  do  exist  for  suitably  limited  domains  of  clutter.  This 
agrees  with  the  near-ideal  performance  that  has  been  observed  in  computer  simulations  using 
TDMA  (or  DDMA)  waveforms  in  scenarios  having  limited  domains  of  clutter.  For  TDMA, 
the  cases  in  which  the  clutter  is  suitably  limited  are  defined  in  terms  of  the  MIMO  radar 
system  design,  involving  the  velocity  of  the  platform,  the  center  frequency,  the  azimuth  and 
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elevation  footprint  of  the  transmit  and  receive  antennas,  the  number  of  transmitters,  and  the 
duty  cycle.  For  clutter  exceeding  such  limited  domains,  the  best  design  of  MIMO  waveforms 
for  minimal  degradation  of  clutter  cancellation  and  other  beneficial  properties  remains  an 
open  issue.  Airborne  experimental  data  will  be  used  to  assess  the  resulting  performance. 
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1.  INTRODUCTION 


Coherent,  multiple-input  multiple-output  (MIMO)  extensions  to  radar  systems  offer  a  number  of 
performance  advantages  over  more  conventional  approaches  [1-1 1].  For  ground  moving  target  indication 
(GMTI),  the  benefits  can  include  lower  minimum  detectable  velocity  (MDV)  [1,  3,  8,  9]  and  more 
accurate  target  direction  estimation  [1,  3].  In  part,  these  benefits  are  afforded  by  the  superior  clutter 
cancellation  offered  by  the  long  filled  virtual  array  that  MIMO  can  yield,  while  using  relatively  few 
physical  antenna  elements  [1,3,  8].  In  addition,  the  use  of  uncorrelated  waveforms  [12-14]  for  coherent 
MIMO  GMTI  diffuses  the  transmitted  energy  uniformly  over  a  broad  azimuthal  coverage  region,  yielding 
a  long  effective  coherent  integration  time  [3,  8].  This  enhances  the  above  benefits  and  also  provides  finer 
Doppler  resolution.  For  applications  in  which  the  geometry  and  MIMO  radar  system  design  result  in 
ground  clutter  occupying  a  suitably  limited  extent  of  range  and  Doppler,  waveform  sets  such  as  time- 
division  multiple  access  (TDMA)  or  Doppler-division  multiple  access  (DDMA)  can  be  used  to  realize  the 
above  benefits  [3,  15].  However,  it  has  been  noted  that  the  clutter  cancellation  is  less  effective  when  the 
clutter  has  a  greater  range-Doppler  extent  [3,  15,  16].  Various  alternative  low-correlation  waveform  sets 
have  failed  to  mitigate  this  problem  [3,  15,  16]. 

In  this  report  it  is  shown  analytically  that  the  extra  difficulty  of  clutter  cancellation  referred  to 
above  is  attributable  to  an  unavoidable  increase  in  the  number  of  MIMO  degrees  of  freedom  required  to 
cancel  the  ground  clutter,  when  the  clutter  occupies  an  extended  domain  of  range  Doppler.  To  explain  this 
effect,  a  pictorial  discussion  is  given  in  this  Introduction,  showing  how  matched  filters  (MF)  for  different 
waveforms  effectively  weight  the  complex  clutter  reflections  from  the  ground  with  different  correlation 
functions  versus  range.  The  different  correlation  functions  result  in  decorrelation  of  the  outputs  of  the 
filters  matched  to  the  different  waveforms.  In  turn,  this  causes  degradation  of  the  ability  to  cancel  the 
clutter  returns  when  the  matched-filter  outputs  are  adaptively  weighted  for  that  purpose.  Time-domain 
computer  simulation  results  are  presented  [16],  showing  how  the  signal-to-interference-plus-noise  ratio 
(SINR)  for  MIMO  GMTI  is  less  than  ideal  for  extended  range  domains  of  ground  clutter.  Specifically, 
practical  waveforms  chosen  to  have  low  cross-correlation  functions  over  extended  range  domains  are  seen 
to  result  in  significantly  worse  SINR  than  do  hypothetical  ideal  waveforms.  The  ideal  waveforms  are 
simulated  as  having  matching  autocorrelation  functions  and  zero  cross-correlation  functions. 

Figure  1.1  depicts  the  antenna  footprint  on  the  ground,  for  an  ordinary  single-transmitter  GMTI 
radar.  The  regions  of  clutter  on  the  ground  experiencing  the  highest  system  gain  will  be  examined.  The 
radar  is  assumed  to  fly  horizontally  across  the  bottom  of  the  diagram.  The  range  and  azimuth  variations 
on  the  ground  are  as  indicated.  The  location  of  a  hypothesized  target  is  shown  in  blue.  The  target  has  a 
hypothesized  range,  azimuth,  and  radial  velocity.  At  the  lower  left,  the  initial  stages  of  the  receive  signal 
processing  are  depicted.  The  signal  processing  is  intended  to  focus  each  antenna  output  on  the  range  and 
Doppler  of  the  hypothesized  target.  This  causes  specific  regions  on  the  ground  to  experience  high  system 
gain.  The  matched  filters  (MFs)  all  focus  on  the  target  range,  and  their  peak  response  is  shown  in  green  in 
the  figure.  The  Doppler  beam  formers  focus  the  MF  outputs  on  the  hypothesized  Doppler  of  the  target. 
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and  their  peak  response  on  the  stationary  ground  is  shown  in  red  (a  tfcray”).  Note  that  the  ground 
experiencing  the  highest  gain  is  not  at  the  same  azimuth  angle  as  the  target.  This  is  because  the  target  is 
hypothesized  to  be  moving  radially,  so  the  ground  in  the  direction  of  the  target  does  not  have  the  same 
Doppler  shift  as  the  target  does;  hence  it  does  not  receive  peak  gain  as  the  target  does.  When  the 
hypothesized  Doppler  of  the  target  is  changed,  the  azimuth  angle  of  the  clutter  isoDoppler  to  the  target 
also  changes,  approaching  the  angle  of  the  target  as  the  hypothesized  radial  velocity  of  the  target  relative 
to  the  ground  approaches  zero.  The  ground  region  receiving  the  highest  overall  gain  is  shown  in  gold;  it  is 
the  intersection  of  the  peak  range  response  of  the  MFs  and  the  peak  Doppler  response  of  the  Doppler 
beam  focused  on  the  target. 
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Figure  1.1.  Primary'  ground  clutter  regions  for  GMTl  (single-transmitter  case). 


The  region  of  ground  clutter  that  is  the  focus  of  the  analysis  is  shown  in  Figure  1.2.  This  is  the 
clutter  ray  (in  red)  isoDoppler  to  the  hypothesized  target.  Here  is  depicted  (in  green)  the  (complex) 
correlation  response  of  the  MFs,  peaking  at  the  hypothesized  target  range  [17,  pp.  32-33].  The  correlation 
response  weights  the  reflectivities  of  patches  of  ground  on  the  clutter  ray.  This  occurs  in  each  antenna 
output  channel,  as  mentioned  above. 
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Figure  1.2.  Ground  clutter  region  for  analysis  (single-transmitter  case). 


The  next  step  in  the  processing  is  to  combine  the  outputs  of  the  antennas  in  an  adaptive  manner,  to 
cancel  the  clutter.  This  is  depicted  in  Figure  1.3.  The  array  gain  pattern  versus  azimuth  angle  is  depicted 
in  red  at  the  bottom  of  the  diagram.  The  gain  is  lower  in  the  direction  of  the  clutter  ray  and  higher  in  the 
direction  ol  the  target.  If  the  hypothesized  target  had  a  lower  velocity  (relative  to  the  ground)  than  shown 
in  Figure  1 .3,  then  the  clutter  isoDoppler  to  the  target  would  be  at  a  closer  angle  to  the  target.  This  would 
reduce  the  ability  ol  the  array  to  form  an  adaptive  null  on  the  angle  of  the  clutter  ray  while  maintaining 
gain  on  the  angle  of  the  target.  Figure  1 .4  shows  computer  simulation  results  [  1 6]  for  the  SINK  that  would 
be  achieved  in  a  sample  scenario  (see  details  on  figure).  The  abscissa  shows  the  normalized  Doppler  [18, 
p.  13]  of  the  hypothesized  target,  relative  to  the  ground.  The  ordinate  shows  the  S1NR  loss,  in  dB.  relative 
to  a  clutter-free  case.  The  loss  is  worst  at  0  normalized  Doppler,  when  the  clutter  isoDoppler  to  the  target 
is  at  the  same  angle  as  the  target  and  thus  cannot  be  nulled  relative  to  the  target  through  adaptive 
weighting  of  the  antenna  outputs.  The  result  shown  is  basically  the  same  for  two  different  radar  system 
structures.  One  is  the  single-input  multiple-output  (SIMO)  structure  with  a  25-element  receive  array  at 
(k/2)  (i.e.,  critical)  spacing.  The  second  is  a  5-transmitter,  5-receiver  MIMO  array,  with  ideal  waveforms 
simulated,  as  noted  on  the  figure.  The  virtual  array  for  MIMO  also  has  25  elements,  at  (X./2)  spacing2.  The 
simulation  results  are  basically  the  same  whether  the  critically  sampled  virtual  array  is  achieved  by  using 


The  locations  of  the  antenna  elements  of  the  MIMO  virtual  array  are  the  result  of  the  convolution  of  the  transmit 
element  locations  with  the  receive  element  locations  [1,3,  8]. 
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receivers  at  (A./2)  spacing  and  transmitters  at  (5X./2)  spacing,  or  the  reverse.  The  reasons  that  these  three 
cases  are  equivalent  will  be  discussed  in  Section  6,  relative  to  (6.28).  Incidentally,  if  the  ground  clutter 
had  an  extended  Doppler  domain  (e.g.,  due  to  a  high-velocity  platform),  then  there  would  be  multiple  rays 
of  clutter  on  the  ground  (i.e.,  Doppler  ambiguities).  Those  rays  would  each  require  cancellation, 
increasing  the  rank  of  the  MIMO  clutter  covariance  matrix.  This  topic  arises  again  in  Section  8.  In 
Figure  1.3,  the  ability  to  cancel  the  clutter  effectively  depends  on  the  antenna  outputs  having  high 
channel-to-channel  correlation  for  the  clutter.  The  fact  that  the  channels  utilized  identical  MFs  in  the 
SIMO  case  has  ensured  that  the  weighting  functions  on  the  ground  are  identical  on  the  clutter  ray.  This 
maintains  high  correlation  between  channels,  hence  effective  clutter  cancellation. 


Channel-to-channel  correlation 
of  clutter  is  high,  because  all 
channels  have  the  same 
weighting  of  the  clutteron 
the  ray 

Clutter  cancellation  benefits 
from  high  channel-to-channel 
correlation 
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Figure  1.3.  Adaptive  gain  pattern  (single-transmitter  case). 
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Parameters 

•  T  =  5,  R  =  5  MIMO  system 

•  Critically  sampled  virtual  array 

•  Ideal  MIMO  waveforms 

•  2  GHz  center  frequency 
•1.5  kHz  PRF 
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•  40  m/s  platform  speed 
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Figure  1.4.  Computer  simulation  results  using  idea l  MIMO  waveforms  in  a  sample  scenario.  Mote  that  all 
simulation  results  in  this  report  utilize  2  pulse-repetition-interval  (PRI)  staggers  /IS,  p.  125],  for  somewhat 
improved  performance  over  simple  post-Doppler,  angle-only  adaptivity,  which  is  analyzed  in  the  text 


In  contrast,  the  situation  tor  MIMO  reception  using  realistic  waveforms  is  now  examined. 
Figure  1.5  depicts  reception  for  a  3-transmitter  case.  More  details  of  the  signal  processing  structure  are 
given  in  Section  3.  Each  receiving  antenna  is  followed  by  three  MFs,  one  for  each  waveform.  For  any 
antenna,  the  weighting  of  the  clutter  at  the  output  of  MF1  depends  on  the  autocorrelation  function  of 
waveform  #1  and  on  its  cross-correlation  functions  with  waveforms  #2  and  #3.  Unless  the  waveforms  are 
the  same,  it  turns  out  (as  proven  in  Section  6)  that  the  clutter  on  the  ray  is  not  weighted  identically  at  the 
three  MF  outputs,  causing  a  loss  of  correlation  (this  can  be  avoided  when  the  clutter  occupies  only  a 
limited  domain  of  range  Doppler,  as  discussed  in  Section  8).  This  limits  the  ability  to  cancel  the  clutter.  In 
fact,  if  the  three  waveforms  were  fully  uncorrelated  (at  all  delays),  then  the  MIMO  clutter  covariance 
would  have  3  equal  eigenvalues,  requiring  3  degrees  of  freedom  to  cancel  the  clutter  (see  Section  4).  This 
is  in  contrast  to  the  I -transmitter  case  discussed  above,  in  which  all  the  MF  outputs  would  have  high 
correlation,  and  the  clutter  covariance  in  theory  could  have  rank  I,  as  discussed  in  Section  3.  The 
performance  impact  of  using  realistic  low-correlation  MIMO  waveforms  is  seen  in  Figure  1.6  1 16]  The 
result  for  a  MIMO  system  with  a  specific  set  of  low-correlation  waveforms  (“Wolf  Whistles  ')  have  been 
added  to  the  results  for  ideal  waveforms  from  Figure  1.4.  The  width  of  the  SINR  notch  around  0 
normalized  Doppler  has  clearly  increased.  The  extra  wiggle  in  the  SINR  loss  curve  at  normalized 
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Doppler  near  0. 1  is  due  to  an  extra  target  in  the  training  data  for  the  adaptive  weights  '.  A  number  of  low- 
correlation  waveform  sets  have  been  tried  [16],  with  similar  results  to  those  in  Figure  1.6.  In  this 
particular  case,  the  five  receivers  had  ( X/2 )  spacing,  and  the  transmitters  had  (5A./2)  spacing.  In  Section  4, 
the  (significant)  impact  of  reversing  the  transmitters  and  receivers  for  realistic  waveforms  will  be 
discussed. 
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Figure  1.5.  Reception  for  Ml  MO  GMT l  (3 -transmitter  case). 


The  time-domain  GMT1  simulation  included  clutter,  receiver  noise,  and  also  5  moving  targets,  having  various 
ranges,  azimuths,  and  velocities.  The  sample  covariance  matrix  used  to  fonn  the  adaptive  weights  was  trained  by 
sampling  the  outputs  of  the  Doppler  beam  formers  (see  Figure  1.5)  at  the  normalized  Doppler  frequency  of  the 
target  hypothesis.  Training  samples  were  taken  at  a  set  of  300  range  bins,  not  including  the  range  bin  of  the  target 
hypothesis  (to  prevent  self-nulling).  In  the  case  shown  in  Figure  1.6,  an  actual  target  was  located  in  one  of  the  range 
bins  used  for  training.  When  the  normalized  Doppler  frequency  of  the  hypothesized  target  equaled  that  of  the  actual 
target,  the  sample  covariance  matrix  contained  significant  power  from  the  actual  target.  This  forced  an  array  null  to 
be  formed  on  the  angle  of  the  actual  target,  in  addition  to  the  null  on  the  clutter  ray  isoDoppler  to  the  target 
hypothesis.  The  formation  of  two  nulls  required  two  degrees  of  freedom,  for  ideal  waveforms.  However,  more 
degrees  of  freedom  were  required  for  realistic,  low-correlation  waveforms,  causing  the  extra  wiggle  in  the  S1NR  loss 
curve  in  Figure  1.6,  at  the  normalized  Doppler  frequency  of  the  actual  target  in  the  training  data. 
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•  2  GHz  center  frequency 
•1.5  kHz  PRF 

•  200  pulse  CPI 

•  40  m/s  platform  speed 

•  2  PRI  staggers 


Figure  1.6  Realistic,  low-correlation  waveforms  (“Wolf  Whistles”)  yield  poorer  S1NR  performance  than  idea! 
MIMO  w  aveforms. 


The  points  made  in  this  introductory  section  are  quantified  analytically  in  Sections  2  through  6. 
Sections  7  and  8  introduce  additional  topics.  Specifically,  the  following  is  a  section-by-section  synopsis 
of  the  discussion  and  results  in  this  report: 

1 .  Introduction — this  section. 

2.  The  assumptions,  limitations,  and  applicability  of  the  ensuing  analytical  results  are  discussed. 
The  analysis  considers  cases  in  which  the  MIMO  virtual  array  is  critically  spaced  and  filled.  In 
particular,  the  results  consider  only  the  clutter  on  the  ray  isoDoppler  to  the  hypothesized  target.1 * * 4 
The  most  familiar  interpretation,  and  that  used  in  this  Introduction,  would  be  that  the  analytical 
results  apply  to  a  pulsed-Doppler  paradigm,  where  each  waveform  consists  of  a  sequence  of 
identical  pulses  over  the  coherent  processing  interval,  and  Doppler  filtering  (to  focus  on  the 
hypothesized  target  Doppler)  is  performed  after  matched  filtering  and  prior  to  adaptive 
weighting  for  clutter  cancellation  (i.e.,  post-Doppler  clutter  cancellation  [18]).  However,  a 
broader  interpretation  is  that  the  waveforms  are  free  to  be  of  any  form  (pulsed  or  not),  and  the 
analysis  here  shows  the  statistical  structure  of  the  clutter  arising  from  just  the  most  critical 
region  on  the  ground  (i.e.,  isoDoppler  to  the  target). 


1  In  some  system  designs  with  low-velocity  platforms,  there  can  occur  hypothesized  target  velocities  tor  which  no 

isoDoppler  clutter  ray  exists  on  the  ground  [18.  p.  27].  Clutter  cancellation  may  not  be  required  in  such  cases  Also, 

if  there  are  multiple  rays  isoDoppler  to  the  target  (i.e.,  Doppler  ambiguities)  within  the  azimuthal  cov  erage  region, 
clutter  cancellation  requires  more  degrees  of  freedom  in  order  to  place  antenna  array  nulls  in  multiple  directions. 

The  analysis  herein  considers  only  cases  in  which  a  single  isoDoppler  clutter  ray  occurs. 
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3.  A  mathematical  model  is  developed  for  the  response  of  the  M1MO  receive  antennas  and  MFs  to 
the  ray  of  clutter  on  the  ground  isoDoppler  to  the  hypothesized  target.  Each  of  the  R  receivers 
is  followed  by  T  matched  filters,  corresponding  to  the  T  transmit  waveforms.  The  RT  outputs, 
sampled  at  the  range  of  the  hypothesized  target  and  Doppler  filtered  to  focus  on  the 
hypothesized  target  velocity,  are  adaptively  weighted  to  cancel  the  clutter  on  the  isoDop  while 
maintaining  gain  on  the  hypothesized  target  azimuth  angle.  The  study  considers  the  rank  and 
eigenstructure  of  the  RT  x  RT  covariance  matrix  between  the  various  received,  processed 
outputs,  as  resulting  just  from  the  clutter  on  the  isoDop  ray  on  the  ground.  When  more 
eigenvalues  are  positive,  more  degrees  of  freedom  potentially  need  to  be  expended  to  cancel  the 
ground  clutter.  Any  performance  degradation  from  training  the  adaptive  system  to  learn  this 
covariance  is  ignored  here  (i.e.,  the  covariance  is  assumed  to  be  known).  As  noted  above,  Unite 
training  data  is  used  in  the  computer  simulation  results.  In  this  section,  it  is  shown  that  for 
uncorrelated  signals,  the  frequency  spectra  cannot  overlap  and  the  autocorrelation  functions  are 
orthogonal  to  one  another,  which  causes  decorrelation  of  the  ground  clutter  returns  at  the 
outputs  of  the  MFs  for  the  different  signals.  It  follows  that  ideal  waveform  sets  do  not  exist 
over  an  extended  domain  of  range.  This  last  fact  was  already  known  [19],  but  the  extent  of  its 
impact  on  the  ability  to  cancel  clutter  was  not  fully  realized.  The  MIMO  response  to  a  point 
target  is  shown  to  have  the  familiar  Vandermonde  virtual-array  characteristic  of  MIMO,  while 
the  response  to  clutter  can  have  multiple  modes.  Also,  it  is  shown  how  single-input  multiple- 
output  (SIMO)  GMTI  yields  a  clutter  covariance  having  only  one  nonzero  eigenvalue. 

4.  Under  simple  assumptions,  it  is  shown  that  the  MIMO  clutter  covariance  has  T  equal,  positive 
eigenvalues  whenever  the  T  waveforms  are  uncorrelated.  This  occurs  independent  of  the 
direction  of  the  ray  of  clutter  isoDoppler  to  the  target.  This  phenomenon  is  a  worst  case  for 
clutter  cancellation.  Also,  it  is  shown  that  in  this  case  the  receive  array  alone  performs  all 
clutter  nulling  for  MIMO,  for  each  transmitter  separately.  This  makes  MIMO  and  SIMO 
equivalent  for  clutter  cancellation,  in  this  respect. 

5.  For  the  case  where  T  =  2,  a  penalty  function  is  derived,  evidencing  quadratic  growth  of  the 
product  of  the  highest  two  eigenvalues  of  the  MIMO  clutter  covariance  whenever  the  two 
signals  are  not  proportional.  This  result  relates  directly  to  the  size  of  the  second  eigenvalue  of 
the  covariance,  which  requires  an  extra  degree  of  freedom  for  clutter  cancellation.  The  penalty 
emphasizes  waveform  differences  at  frequencies  where  the  energy  density  spectrum  of  either 
waveform  is  high,  suggesting  the  use  of  flat-spectrum  waveforms.  The  following  statement  is 
also  proven:  For  T  =  2,  the  MIMO  clutter  covariance  is  of  rank  1  or  0  for  all  (|)  if  and  only  if 
(iff)  the  two  waveforms  are  proportional.  Here,  <))  is  defined  as  the  phase  difference  between 
adjacent  transmitters,  as  would  be  seen  on  the  ray  of  clutter  isoDoppler  to  the  target. 

6.  For  arbitrary  T,  it  is  shown  that  the  MIMO  clutter  covariance  is  of  rank  1  or  0  for  all  (j)  iff  the 
T  x  T  waveform  cross-correlation  matrix  as  a  function  of  delay  satisfies  the  equality 
C(r)  =  C(0)/(r),  for  all  relevant  values  of  x,  for  some  scalar  (autocorrelation)  function, 


8 


/(r) .  This  property  is  called  delay  similarity  of  the  waveform  cross-correlation  matrix.  Using 
the  above  equivalence  proves  that  the  statement  following  the  colon  in  (5)  above  is  true  for  any 
value  of  T.  Further,  the  properties  of  analytic  functions  then  are  used  to  extend  the  statement 
finally  to  the  following  pair  of  statements:  (1)  For  any  T,  if  all  the  waveforms  are  proportional, 
then  the  MIMO  clutter  covariance  is  of  rank  1  or  0  for  all  <(),  and  (2)  for  any  T,  if  any  pair  of  the 
T  waveforms  are  not  proportional,  then  the  MIMO  clutter  covariance  can  have  rank  I  or  0  only 
for  isolated  values  of  <j). 

7.  Three  different  potential  techniques  are  examined  briefly  for  improving  the  cancellation  of 
clutter  over  extended  domains.  The  three  techniques  are  adaptively  weighted  slow-time  taps 
(i.e.,  Space-Time  Adaptive  Processing  fSTAP]  [18]),  adaptively  weighted  fast-time  taps  and 
nonconventional  reception  (i.e.,  reception  not  constrained  to  use  matched  filters  for  the  T 
waveforms).  For  uncorrelated  MIMO  waveforms,  it  is  shown  that  none  of  these  techniques 
addresses  the  basic  decorrelation  problem  of  the  MIMO  clutter  covariance.  Thus,  while  these 
techniques  could  improve  performance  somewhat,  they  are  not  expected  to  restore  near-ideal 
performance. 

8.  For  clutter  occupying  only  a  suitably  limited  domain  of  range  Doppler,  it  is  already  known  that 
near-ideal  clutter  cancellation  can  be  achieved  through  appropriate  waveform  design,  as  noted 
above.  Using  the  example  of  TDM  A  waveforms,  one  can  show  here  that  it  is  possible  to 
achieve  the  delay-similarity  property  defined  in  (6)  above  when  the  clutter  occupies  only  a 
suitably  limited  domain,  which  implies  that  the  MIMO  clutter  covariance  will  have  rank  1  or  0 
for  all  (j).  Thus,  ideal  sets  of  waveforms  do  exist  for  suitably  limited  domains  of  clutter.  This 
agrees  with  the  near-ideal  performance  that  has  been  observed  in  computer  simulations  using 
TDM  A  (or  DDMA)  waveforms  in  scenarios  having  limited  domains  of  clutter  [16].  For 
TDMA,  the  cases  in  which  the  clutter  is  suitably  limited  are  defined  in  terms  of  the  MIMO 
radar  system  design,  involving  the  velocity  of  the  platform,  the  center  frequency  ,  the  azimuth 
and  elevation  footprint  of  the  transmit  and  receive  antennas,  the  number  of  transmitters,  and  the 
duty  cycle.  For  clutter  exceeding  such  limited  domains,  the  best  design  of  MIMO  waveforms 
for  minimal  degradation  of  clutter  cancellation  and  other  beneficial  properties  remains  an  open 
issue.  Airborne  experimental  data  will  be  used  to  assess  the  resulting  performance. 

9.  Summary. 

Note  that  [15]  examined  the  difficulties  of  waveform  selection  posed  by  clutter  of  unlimited  range 
or  Doppler  extents  for  a  restricted  situation  (i.e.,  2  transmitters,  1  receiver)  and  suggested  some  MIMO 
waveforms  for  use  in  applications  in  which  the  range  and  Doppler  extents  of  the  clutter  are  limited.  The 
work  herein  extends  the  conclusions  of  [15]  to  an  arbitrary  number  of  transmitters  and  receivers,  and 
quantifies  the  diffusion  of  power  into  multiple  eigenvalues  of  the  MIMO  clutter  covariance. 
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2.  ASSUMPTIONS,  LIMITATIONS,  AND  APPLICABILITY  OF  ANALYTICAL 

RESULTS 


To  generate  insight  herein,  the  GMT1  scenario  and  the  M1MO  design  are  kept  simple,  and  only  the 
primary  contributing  clutter  is  considered  in  the  analysis.  The  MIMO  waveforms  are  considered  to  he 
narrowband,  and  only  far-field  ground  reflections  are  considered.  The  transmitters  are  assumed  to 
transmit  for  a  single  coherent  processing  interval.  The  platform  is  assumed  to  fly  in  a  direction  parallel  to 
the  linear  array  of  transmitters  and  receivers  (sensors)  with  no  yaw  or  pitch  rotation.  This  keeps  the 
geometry  of  the  clutter  simple  on  the  ground,  in  terms  of  its  Doppler  shift  versus  the  cone  angle  off  the 
sensor  array  axis.  The  clutter  having  the  same  Doppler  shift  as  the  target  hypothesis  is  approximated  as 
forming  a  thin  locus  on  the  ground,  at  constant  cone  angle  (it  will  be  called  a  ray).  Only  that  clutter  ray  is 
considered  in  the  analysis.  This  same  approximation  was  made  in  [15];  it  allows  attention  to  be  focused 
on  the  limitations  to  clutter  cancellation  arising  merely  from  the  range  extent  of  the  clutter  on  the  primary 
ray.  As  pointed  out  in  the  Introduction,  the  analysis  considers  only  the  cases  in  which  there  is  exactly  one 
such  ray  of  clutter  isoDoppler  to  the  target,  for  the  target  velocity  of  interest.  To  obtain  insight,  one 
approximates  the  clutter  ray  as  extending  on  both  sides  of  the  hypothesized  target  range,  to  distances 
beyond  that  corresponding  to  the  finite  coherent  processing  interval  of  the  radar.  For  the  mathematical 
model,  this  is  equivalent  to  an  unbounded  extent  of  clutter  in  range,  both  positive  and  negative,  from  the 
range  hypothesized  for  the  target.  In  Section  8,  however,  clutter  of  bounded  extent  is  discussed,  as  could 
occur  because  of  the  impact  of  antenna  footprints  or  of  the  finite  propagation  horizon  of  the  earth.  As 
mentioned  in  Section  I,  near-ideal  SINR  performance  can  be  achieved  when  the  clutter  is  of  sufficiently 
bounded  extent.  Section  8  will  show  why  this  is  true,  in  terms  of  the  conceptualization  developed  herein. 

The  physical  array  of  sensors  has  equally  spaced  transmitters  and  equally  spaced  receivers, 
configured  to  yield  a  critically  sampled  virtual  array  (see  Section  3).  F.ilher  the  T  transmitters  or  R 
receivers  can  be  placed  farther  apart  (details  will  be  given  below).  For  the  MIMO  system,  each  receiver  is 
followed  by  a  matched  filter  for  each  transmitted  waveform,  sampled  at  the  range  bin  of  the  target 
hypothesis.  This  yields  RT  outputs  for  the  range  bin. 

Because  attention  is  restricted  to  a  thin  ray  of  clutter  isoDoppler  to  the  target,  the  analysis  applies 
when  viewed  from  either  of  two  perspectives.  On  the  one  hand,  it  can  be  assumed  that  a  repetitive  pulsed- 
Doppler  radar  waveform  structure  has  been  imposed  on  all  the  waveforms,  and  that  the  matched-filter 
outputs  have  each  been  Doppler  filtered  to  leave  just  the  clutter  isoDoppler  to  the  target  Angle-only 
clutter  cancellation  is  then  applied  to  the  RT  outputs  for  the  Doppler  bin.  The  analytical  results  herein 
then  apply  to  the  (RT  x  RT)  MIMO  clutter  covariance.  The  possibility  for  performance  improvement  in 
practice  by  using  Space-Time  Array  Processing  (STAR)  [18]  is  discussed  in  Section  7.  On  the  other  hand, 
the  restriction  of  consideration  to  just  the  thin  ray  of  clutter  isoDoppler  to  the  target  can  be  seen  as 
providing  free  Doppler  filtering,  hence  allowing  a  broader  class  of  waveforms  to  be  considered.  Again, 
the  clutter  cancellation  is  based  on  the  RT  outputs,  and  the  analytical  results  herein  apply  to  the 


covariance  matrix  of  those  outputs.  Seen  from  this  broader  perspective,  the  negative  predictions  stated  in 
Section  1  are  even  more  impressive. 
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3.  MATHEMATICAL  MODEL  FOR  MIMO  GMTI  OBSERV  ATIONS 


In  this  section  is  developed  a  mathematical  model  related  to  that  utilized  in  [15],  but  modified  to 
highlight  the  issue  of  the  clutter  rank.  The  simple  linear  antenna  array  layout  shown  in  Figure  3.1  is 
utilized.  The  T  transmitters  are  equally  spaced  by  (A72).  The  R  receivers  are  spaced  by  (T  X/2),  resulting  in 
a  critical  QJ2)  spacing  for  the  virtual  array  of  RT  elements.  The  leftmost  antenna  is  used  as  a 
transmit/receive  element  and  is  the  phase  center  used  for  the  mathematical  modeling.  The  platform  is 
assumed  to  travel  in  a  direction  parallel  with  the  array  axis,  for  simplicity.  Figure  3.1  illustrates  the  case 
of  (T  -  3,  R  4),  but  the  analysis  is  performed  for  arbitrary  values  of  T  and  R.  and  holds  equally  if  the 
receivers  were  made  closer  together  and  the  transmitters  farther  apart.  To  focus  on  the  clutter  rank  issue, 
the  analysis  is  performed  by  including  just  the  clutter  on  the  ray  that  is  isoDoppler  to  the  target  (see 
Section  2).  As  noted  in  the  Introduction,  the  analysis  treats  only  the  cases  in  which  there  is  exactly  one 
such  ray.  If  that  clutter  is  at  off-broadside  angle  0  (the  off-broadside  angle  is  the  complement  of  the  cone 
angle),  then  any  point  on  the  clutter  ray  will  observe  (in  the  far  field)  a  phase  between  adjacent 
transmitters  of 


=  jsin#//i  =  ;rsin#,  (3.1) 

due  to  the  k/2  spacing  of  transmitters.  Because  high-speed  targets  may  not  have  isoDoppler  clutter  on  the 
ground,  the  analytical  results  herein  are  simpler  to  state  in  terms  of  4>  rather  than  the  target  velocity. 


Figure  3. 1.  MIMO  transmission  and  reception. 
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The  transmitted  waveforms  are  designated  as  {S|(t),...,sT(t)[,  in  analytical  signal  notation  [20,  pp. 
57-601.  If  the  hypothesized  target  is  at  (slant)  range  r0  from  the  platform,  the  waveforms  would  superpose 
at  that  range  on  the  clutter  ray  proportional  to  the  scalar  time  function 

«-2> 

1=1 

where  the  phasors  represent  the  transmitter  phase  offsets  on  the  clutter  ray5.  For  a  clutter  patch  on  the  ray 
displaced  from  the  target  range  by  r,  the  superposition  in  (3.2)  would  become  instead 


±>,(t-h  +  r\lcyn».  (33) 

/=l 

If  this  patch  has  complex  reflectivity  q(r)  (including  2-way  r-dependent  propagation  loss  and  antenna 
gain),  then  at  the  phase-reference  antenna  (i.e.,  antenna  #1)  the  reception  would  be  proportional  to 

2^,  +  r]/c)e/(M)^(r),  (3.4) 

i=\ 

because  of  this  patch  of  clutter.  The  superposition  of  all  the  clutter  patches  on  the  ray  would  be  observed 
as 


Z  (*  - 2  h + d 7  c)eAi  °Vr)  ■  (15) 

r  1=1 

As  discussed  in  Section  2,  the  focus  here  is  on  the  case  of  an  extended  range  domain  of  clutter,  where  the 
clutter  ray  extends  on  both  sides  of  the  hypothesized  target  range,  to  a  distance  beyond  that  corresponding 
to  the  finite  coherent  processing  interval  of  the  radar.  This  makes  the  sum  over  (r)  in  (3.5)  effectively 
cover  an  infinite  extent  of  range,  both  positively  and  negatively.  The  matched  filter  for  the  kth  waveform, 
at  the  output  of  this  antenna  is  now  considered.  When  that  matched  filter  is  sampled  at  a  time  of  (2ro/c), 


The  expression  in  (3.2)  is  valid  when  the  waveforms  {s,(t)(  represent  the  various  transmissions  over  the  full 
coherent  processing  interval.  For  repetitive,  pulsed-Doppler  radar,  it  is  convenient  to  think  of  the  {Sj(t)}  as  being  just 
the  pulse  waveforms.  For  that  interpretation,  (3.2)  should  include  range  ambiguities.  However,  in  the  assumed 
geometry  the  range  ambiguities  for  a  given  Doppler  frequency  are  at  the  same  angle-of-arrival.  Thus,  while  the 
ambiguities  contribute  added  clutter  power,  they  do  not  change  the  functional  form  of  the  M1MO  clutter  covariance. 
Thus,  (3.2)  can  serve  for  both  cases. 


14 


corresponding  to  the  hypothesized  target  range,  it  effectively  correlates  its  input  function  with  sk(t  -  2rn/c) 
[1 7.  p.  153].  For  the  clutter  input  in  (3.5),  the  output  of  the  matched  filter  would  be 


.  rjc)  S i>,(/-2[-i  +  r]/c)e *q(r) 


r  i=\ 


Jl. 


r  /=!  J 

=  c„(-2/-/c)e'"-'’'9(r). 


r  /=! 


where  ckl()  is  the  cross-correlation  between  the  k  and  i,h  waveforms,  defined  by 


Cu(r)  =  \s'k(t)st{t  +  Y)dt. 


Equation  (3.6)  can  be  rewritten  as 


=  X  [C*I  (-2  /-Ac),...  ,ckT  (-2r/e)]vr  (<p)q  (r), 
where  vT((|))  is  the  Vandermonde  vector  of  the  transmitter  phase  offsets  on  the  clutter  ray: 

vr(^)  =  [ . 

Thus,  coming  from  the  T  matched  (liters  following  antenna  #1,  the  following  vector  is  obtained: 


T  A 

-i 


V~r  J 


w  here  C(  )  is  the  matrix  of  waveform  cross-correlations: 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


[c(r)L,=c*<(r)- 


(3.11) 
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The  receivers  arc  spaced  apart  by  (TA72);  hence  the  outputs  of  the  T  matched  filters  following  the  kIh 
receiver  are  the  same  as  the  vector  (3.10),  multiplied  by  the  scalar  phasor  exp[j(k  -  l)T<j)].  When  the 
T-vector  outputs  of  the  R  receivers  are  stacked  as  shown  in  Figure  3.1,  one  obtains  the  RT-vector — the 
R  -  4,  T  =  3  case  is  examined  first,  to  make  a  point  after  (3.14): 


1  0 

0  1 

0  0 

eiM  0 

0  e)M 


0 

0 

1 

0 

0 

eiM 


£C(-2r/c)vr(<%(r). 


eJ9*  0  0 

0  em  0 

0  0  e m 


(3.12) 


For  general  R  and  T,  (3.12)  becomes 


2  =  ®  1T\  Yj C(~2r  1  c)Y.M)<l(r\ 


(3.13) 


where 


1 


ei(R-\W 


(3.14) 


and  ®  denotes  a  Kronecker  product  [21,  p.  21],  It  is  observed  that  the  (12  x  3)  matrix  on  the  right-hand 
side  (RMS)  of  (3.12)  has  three  orthogonal  columns,  for  T  -  3.  This  wants  of  the  possibility  of  a  rank-T 
covariance  for  z,  as  will  be  encountered  below,  when  the  waveforms  are  uncorrelated. 
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As  an  important  aside,  it  is  useful  to  use  (3.13)  to  derive  the  MIMO  response  to  a  hypothesized 
point  target.  The  target  would  be  at  r  -  0.  For  equal-power  waveforms  uncorrelatcd  at  0  delay,  C(0) 
becomes  Cn(0)  IT.  Then  (3.13)  becomes,  for  the  target. 


2,  =[v„(«l)®/r]cn(0)&(ft)«. 

<3IS> 

where  <pt  is  the  phase  between  transmitters  at  the  target,  and  q,  is  the  complex  reflectivity  of  the  target. 
The  Vandermonde  vector  in  (3.15)  is  the  convolution  of  the  receiver  and  transmitter  locations,  which  is 
referred  to  as  the  MIMO  virtual  array  [1,3,  8],  Other  values  for  C(0)  would  give  other  responses  in 
(3.15),  but  always  a  single  steering  vector  for  the  target.  The  fact  that  the  clutter  is  extended  in  range 
forces  the  use  of  (3. 13),  with  its  multiple  steering  vectors,  in  place  of  (3. 15)  fora  point  source.  The  impact 
will  be  seen  below. 

The  clutter  patches  are  typically  modeled  as  having  independent  reflectivities.  As  the  patches 
shrink,  (3.13)  is  replaced  by  the  integral 


z  =  [ v*  ((*)  ®  1T  ]  f C (~2r  / c )vr  (4)  <1  ( r )  dr, 

where  q(r)  is  a  complex,  circular  Gaussian  white-noise  random  process  [22,  pp.  65-72],  corresponding  to 
the  “sandpaper  earth”  approximation  for  clutter,  which  is  held  to  be  valid  when  the  granularity  of  the 
clutter  is  coarser  than  the  wavelength  of  the  radar.  Changing  variables  in  (3. 16),  one  obtains 

2:  =  (^)<8>/r  ]  Jc(r)  vr  (^)^(-c-r  /  2)z/r, 

and  because  q(-cr/ 2)  is  still  a  white  noise  (having  a  scaled  power  spectral  density),  it  is  renamed 
u(x)  for  ease: 


^  =  [v«  (<z>)<S>/r]  Jc?(r)vr  (^)w(r)c/r.  (3.17) 

Because  the  columns  of  the  first  matrix  on  the  RMS  of  (3.17)  are  orthogonal,  it  will  be  observed  that  it's 
the  variation  of  the  T  *  T  matrix  C(x)  over  x  that  is  critical.  This  variation  is  constrained  by  the  fact  that 
C(x)  is  the  cross-correlation  between  actual  waveforms,  resulting  in  the  findings  that  were  outlined  in 
Section  1. 
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Continuing  the  modeling,  the  correlation  function  of  u(x)  is  written  as 


£[w(r)w*(r,)]  =  /?(r)^(r-r1),  (3.18) 

where  £>(•)  is  the  Dirac  delta  function.  The  function  p(r)  includes  the  impact  of  variability  of  antenna 
gains,  terrain  cover,  and  propagation  losses.  For  simplicity,  p(t)  will  be  set  to  unity  in  this  analysis.  To 
simplify  the  notation,  one  can  define  the  following: 


(3.19) 


It  is  observed  that  )J  is  an  orthonormal  set  of  T  vectors,  each  having  dimension  RT.  The 

covariance  of  z  from  clutter  on  the  isoDop  ray  has  been  referred  to  as  the  MIMO  virtual  covariance 
matrix  associated  with  the  target  Doppler  bin  [3].  For  brevity,  it  is  referred  to  here  as  the  MIMO  clutter 
covariance.  From  (3.17-19),  it  can  be  written  as 


2(^)= E [zz" 


5(0)[jjc(r)v7(0)^(r-r,)v;/ 

i?(0)//(0)5"(0), 


(r,)dr</r,  BH  (0) 


(3.20) 


where  the  sifting  property  of  the  Dirac  delta  function  has  been  used,  and  the  resulting  integral  has  been 
abbreviated  by 


where 


(3.21) 


G(r,</>)  =  C(r)vT(</>)v"  (r). 


(3.22) 


Because  G[r,</>)  is  (T  *  T),  positive-semidefinite  Hermitian  for  all  (r,0),  H  (0)  is  also  (T  xT), 
positive-semidcfinitc  Hermitian.  Thus,  one  can  write  an  eigendecomposition  for  //  (0)  as 


//(0)  =  F(0)A(0)Fw(0), 


(3.23) 
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where  A(^)  is  diagonal,  non-negative  real.  Because  j^#(^)|l/V/? 
can  use  (3.23)  to  write  the  eigendecomposition  of  Z(y)  as 


has  orthonormal  columns,  one 


2(^)=[5(^)F(4(l/V^)][A(^)^][fi(^)F(^)(l/V^) 


(3.24) 


where  it  is  clear  that  the  columns  of  the  first  factor  in  brackets  on  the  RHS  are  orthonormal.  Thus,  Z(^) 
has  (RT  -  T)  eigenvalues  of  zero,  and  T  eigenvalues  scaled  up  from  those  of  //  [({>)  (some  of  which  may 
also  be  zero): 


R-MhW) 

0 


:(i  =  r+ 1 . RT) 


(3.25) 


The  eigenvalues  are  considered  to  be  written  in  descending  order. 

Some  mathematical  preliminaries  are  now  covered.  First,  as  is  well  known,  the  cross-correlation 
functions  satisfy 


=  J[J st  J  [  j 5,  (r)eM 'dy]dt 

=  | ]s;(f)s<(r)e'2”'s(f-r)Jrdl 

=  \s;U)Sl(/)e'-’"df,  (3.26) 


where  Sk  (  /  )  represents  the  Fourier  transform  of  sk[t).  Thus,  Sk  (f)S,  (/)  is  the  Fourier  transform 
orc„(r): 


f[c„(d]  =  S;(/)S,(/). 


Because  of  the  expressions  in  (3.21-22),  there  is  particular  interest  in  integrals  of  the  form 
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Jh  (rK,  (r)]<fr  =  J[J  5;  (f)s,  (/)€*'•  df\\ S ;  (r)s,  (y)e^dr]  dr 

“if  si(f)si(f)st(y)s'r(r)\je',"'',''dr\drdf 

=|Js;(/)s)(/)s1(r)s;W<5(/-r)^/ 

=|s;(/)s,(/)s,(/)5;(/)d/  <328> 

4*„.  <329> 

Due  to  (3.29),  there  are  many  symmetries  between  the  indices  of  <S>jjkf ,  for  example, 

V  =4%>=®;,r  <3-30> 


Integrals  of  the  form  (3.29)  are  inviting  for  the  application  of  Schwarz’s  inequality,  which  will  be  utilized 
in  Section  5.  In  the  past,  the  Schwartz  inequality  has  been  used  on  quantities  like  those  in  (3.28),  to 
develop  various  bounds  related  to  correlation  properties  of  discrete  sequences  [23,  24]. 

It  is  assumed  herein  that 


c„(0)>0;(^  =  l,...,r).  (3.31) 

Combining  this  with  (3.29)  yields 

=  J|  c„  (r)  |2  dr  >  0;(*  =  1,...,  T).  (3-32) 

It  previously  had  been  pointed  out  by  Keith  Forsythe  of  Lincoln  Laboratory  that  ideal  waveform 
sets  cannot  in  fact  exist  over  extended  domains  of  range,  where  an  ideal  waveform  set  would  consist  of 
waveforms  that  are  pairwise  uncorrelated  for  all  delays  but  have  identical  autocorrelation  functions  [19]. 
Using  the  tools  above,  one  easily  can  see  why  such  sets  cannot  exist.  If  sk  (/)  and  sf  (/)  are  uncorrelated 
at  all  delays,  then 


Jk*c  (r)|2  dT  =  0. 


(3.33) 
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But  then,  by  (3.29), 


<*>«*/  =0- 


(3.34) 


In  particular,  from  (3.29), 


(3.35) 


(3.36) 


So  Sk  (/)  and  S,  (  / )  share  no  frequeney  support.  Now.  one  can  examine  the  mismatch  between  the 
autocorrelation  functions  ckk  ( r)  and  cn  (r) .  To  do  this,  one  can  use  the  generalized  cosine  between  two 
functions  of  i: 


,, 


=  0,  (3.37) 

where  the  last  step  follows  from  (3.36).  Thus,  ckk  ( r )  and  c,Y(r)  are  orthogonal  (hence  maximally 
mismatched).  From  the  discussion  in  Section  I  of  the  matched-filter  weightings  of  the  ground  clutter  on 
the  isoDopppler  ray,  the  drastic  mismatch  of  autocorrelation  functions  for  the  different  waveforms  of  the 
set,  combined  with  zero  cross-correlation,  is  expected  to  cause  decorrelations  in  the  clutter  covariance 
matrix.  That  will  be  quantified  in  Section  4. 

It  is  noted  here  that  if  all  the  waveforms  are  proportional  (i.c..  SIMO),  then  (3.7,  1 1 )  imply  that 

.x,  (/)  =  </,*(/);(,•  =  i....,r), 

C(r)  =  Jc(r)Jw,  forde  Cr,  (3 
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where  c(r)  is  the  autocorrelation  function  of  s(t).  Then  in  (3.22) 


G(z,<p)  =  d  c(r)dH  vT  (^)  ( </>)dc*{r)d " 

=  d\c(T)d"vT(</>)\2  dH , 


(3.39) 


which  has  at  most  one  nonzero  eigenvalue,  whose  only  potential  eigenvector  (r/)  is  constant  in  direction 
over  the  domain  of  x.  Thus,  H  [</>),  in  (3.21)  has  rank  0  or  1.  In  fact,  if  waveforms  could  be  designed  to 
be  dependent  on  the  hypothesized  value  of  <j),  it  would  be  possible  to  make  //(^)  have  rank  0 — by 
choosing  d  in  (3.39)  to  be  perpendicular  to  vT  (</>) — making  Z(^)  also  be  of  rank  0,  by  (3.25).  This  is 
just  transmit  nulling  of  the  clutter  ray.  In  this  work,  it  is  assumed  that  a  single  set  of  waveforms  is 
employed  for  all  hypothesized  values  of  4>,  which  rules  out  that  possibility.  Note  that  (3.38)  is  just  a 
special  case  of  delay  similarity  of  the  waveform  cross-correlation  matrix  (Condition  C  of  Figure  6.1), 
which  always  guarantees  that  H  {  (fi)  has  rank  1  or  0. 

Returning  to  the  mathematical  preliminaries,  from  (3.9,  1 1,  22)  one  obtains 


G(r,^)  = 


Zci 


1=1 


Zcr,(0 


,/<*('- 0 


_  i=I 


i=l 


~0 


1=1 


(3.40) 


[G(r^)L=Z  2>*/(rK*(0e 


M/-*) 


*=i  /=i 


(3.41) 


Hence,  from  (3.21, 29) 


["WL-Z  I**"*1® 


mwk  * 


k=\  /=! 


(3.42) 


Equation  (3.42)  will  be  used  in  the  next  section. 
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4.  MEMO  CLUTTER  COVARIANCE  FOR  UNCORRELATED  WAVEFORMS 


When  the  waveforms  arc  pairwise  uncorrelated,  equations  (3.33-36)  hold  for  any  pair  of 
waveforms  from  the  set.  Thus,  for  (3.42), 

=  0,  unless  m  =  i  =  n  =  k.  (4.1) 


But  then,  from  (3.42). 


;  ( m  =  n ) 

;  otherwise 


If  one  merely  assumes  that  in  (3.32) 


O  —  OifVw), 

mmmm  ’  \  /  ’ 


then 


//(^)  =  :rO/7;(V^), 


(4.2) 


(4.3) 


(4.4) 


having  T  equal  eigenvalues.  However,  recall  from  (3.25)  that  the  nonzero  eigenvalues  of  the 

M1MO  clutter  covariance,  are  just  scaled  up  from  those  of  H  [(p) .  Thus,  X(^)  also  has  T  equal,  positive 
eigenvalues,  a  poor  situation  for  clutter  cancellation.  Note  that  this  is  true  for  all  values  of  <|>. 

The  result  in  (4.4)  also  allows  one  to  see  why  the  SINR  loss  curve  as  a  function  of  velocity  has 
different  shapes,  depending  on  whether  the  transmitters  or  receivers  are  spaced  far  apart  physically.  For 
uncorrelatcd  waveforms,  (4.4)  says  that  H  ((/>)  is  a  sealed  identity  matrix.  But  any  set  of  weights,  w,  that 
fully  cancels  the  clutter  must  satisfy 


E 


=  0 


n"X(^)vy  =  0. 


So.  using  (3.20)  and  (4.4).  the  weights  must  satisfy 


(4.5) 
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(4.6) 


w" B(</>)T<S>ItBh  (0)w  =  O. 

T<t>\\B"  (0)w\\2=O. 

Thus,  w  must  be  orthogonal  to  each  column  of  B (</>),  as  exemplified  by  the  first  factor  on  the  RHS  of 
(3.12),  and  written  out  in  (3.19).  However,  these  columns  have  no  common  elements  between 
transmitters,  relying  instead  on  just  the  receiver  outputs  for  a  single  transmitter  to  cancel  one  another. 
Thus,  when  the  receivers  are  close  together,  the  SINR  notch  is  wide  due  to  poor  angular  resolution,  as 
shown  in  Figure  1.6.  When  the  receivers  are  spaced  at  greater-than-critical  spacing,  there  are  multiple  thin 
notches  in  the  SINR  loss  curve  because  of  ambiguities  in  the  receiver  array  response.  The  computer 
simulation  results  for  this  case  are  shown  in  Figure  4.16.  The  fact  that  the  receive  array  alone  is 
responsible  for  clutter  cancellation  when  the  waveforms  are  uncorrelated  makes  MIMO  perform  like 
SIMO  in  this  respect. 


Notation:  (T  spacing,  R  spacing)  in  XI 2  units 


Normalized  Target  Doppler 


Parameters 

•  T  =  5,  R  =  5  MIMO  system 

•  "Wolf  Whistle”  waveforms 

•  2  GHz  center  frequency 
•1.5  kHz  PRF 

•  200  pulse  CPI 

•  40  m/s  platform  speed 

•  2  PRI  staggers 


Figure  4. 1 .  The  blue  curve  repeats  results  from  Figure  1.6,  with  closely  spaced  receivers .  The  red  cur\>e  is  for 
widely  spaced  receivers,  causing  array  ambiguities . 


Recall  from  the  discussion  of  Figure  1 .6  that  there  was  an  extra  target  in  the  training  data,  causing  extra  SINR  loss 
at  a  hypothesized  normalized  target  Doppler  of  about  0.1.  That  extra  loss  affects  both  array  layouts  in  Figure  4.1, 
due  to  the  use  of  realistic,  low-correlation  waveforms,  which  require  extra  degrees  of  freedom  to  attempt  to  cancel 
the  isoDop  clutter  ray  together  with  the  direction  of  the  true  target  in  the  training  data. 
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5.  IMPACT  OF  UNEQUAL  WAVEFORMS  ON  THE  SECOND  EIGENVALUE 
OF  THE  MIMO  CLUTTER  COVARIANCE,  WHEN  T  =  2 


In  this  section,  it  shall  be  seen  that  when  T  =  2  there  is  a  quadratic  penalty  on  the  product  of  the 
highest  two  eigenvalues  of  the  MIMO  clutter  covariance  when  the  two  waveforms  are  not  proportional. 
This  yields  information  about  the  behavior  of  the  second  eigenvalue. 

One  starts  by  writing  out  (3.22),  using  (3.9,  1 1 ); 


G’(r,^)  = 


G  i  (  0  ci  i  ( t )  +  c,  |C|2^ 

j<h  * 

CIlC2l  +C1 

+c12e,V*1  +  cl2c,*2(r) 

^21 

<"2l<a  1  +  C2\C\2e 

^*21^*21  ^  2 

+c22el‘pc'n  +  c22c\2  ( r ) 

+c22e"*'c2\  - 

j4> 


(5.1) 


where  the  argument  (r)  has  been  left  off  of  most  of  the  cy  (r) ,  for  brevity  .  As  indicated  in  (3.21),  one 
must  integrate  G(r,<p)  over  r  to  obtain  H {</)).  Using  (3.29),  one  obtains: 


H(+)  = 


+<t>1211e'*+012I2 

<*>2111  +  <*>2112^ 
+(I)22ll£’t  +  (t>2212 


+(f)122|CiC'  +Ct>|223 

(U  2,21  +  <&H22e  4 

<*>2221*" +<*>2222 


(5.2) 


Recall  from  Section  3  that  G(t,<p)  and  U  {</>)  are  both  positive-semidelimte  Hermitian  matrices.  Thus, 
the  second  eigenvalue  of  H  ((/>)  is  never  negative,  and  is  positive  iff  the  determinant  of  H  (<p)  is 
positive.  The  mathematical  expression  for  the  determinant  is  simpler  than  that  for  the  second  eigenvalue, 
so  only  the  determinant  is  studied  here.  Because  many  of  the  exponentials  in  the  determinant  integrate  to 
0,  one  examines  the  following  expression,  instead  of  studying  various  values  of  <p  individually: 


A  verage  J  det  [  H  ( (p )  J  j  . 

<!> 


(5.3) 


When  the  quantity  in  (5.3)  is  positive,  then  Amin[f/(^) ]  cannot  be  0  for  all  and  so  /,[  X(^)  j  in 
(3.25)  cannot  be  0  for  all  <p ;  thus,  Z(^)  (the  MIMO  clutter  covariance)  cannot  have  rank  I  or  0  for  all 
<p  .  A  quadratic  penalty  on  (5.3)  will  be  derived,  which  applies  when  the  waveforms  are  not  proportional. 
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Continuing  the  analysis,  one  recalls  that 


=  W||W22  -  w,2w2|. 


(5.4) 


If  one  applies  this  to  (5.2)  and  averages  over  (f) ,  (5.3)  becomes 


(5.5) 


To  derive  this  expression,  some  simple  symmetries  similar  to  (3.30)  have  been  used.  The  Schwarz 


inequality  will  be  employed  (multiple  times)  on  the  RHS  of  (5.5).  The  version  employed  includes  a 
remainder  term,  and  is  equivalent  to  the  Pythagorean  theorem,  which  is  now  discussed. 

The  Pythagorean  theorem  is  now  examined,  as  applied  to  the  L2  space  of  square-integrable, 
complex  functions  (see  [25,  pp.  49-52]).  Its  relation  to  the  Schwarz  inequality  will  also  be  brought  out.  In 
Figure  5.1,  a(f)  and  b(f)  are  complex  functions  of  the  continuous  variable  (0-  In  the  analysis  herein,  (1) 
represents  frequency.  Different  values  of  (0  can  be  considered  as  the  coordinates,  with  a(f)  and  b(f)  being 
infinite-dimensional  vectors  on  these  coordinates,  in  the  vector  space  of  L2.  The  dot  product  in  this  vector 
space  is  analogous  to  that  of  finite-dimensional  vectors: 


(5.6) 


and 


(5.7) 


with  equality  iff  b(  )  is  the  0  function  (sets  of  measure  0  can  be  ignored,  because  it  would  take 
waveforms  of  infinite  duration  to  generate  a  Dirac  delta  function  in  the  frequency  domain).  Define,  when 
(b)  is  of  nonzero  norm. 
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(5.8) 


±M/)M 

II  b  ||2 

This  quantity  is  called  the  projection  of  (a)  onto  (b),  because  Pba  is  parallel  to  (b) — obv  ious  from  (5.8) — 
and  the  residual  error 


e(f)=a(f)-Pha 


(5.9) 


is  perpendicular  to  (b)  (and  hence  is  also  perpendicular  to  Pha): 


[*'(/)>/ 

H 

=  0. 

Thus.  Figure  5. 1  shows  a  right  angle  between  (e)  and  (b).  From  (5.9),  one  may  write  a(f)  as 

a(f)  =  e{f)  +  Pha, 


(5.10) 


(5.11) 


where  the  two  terms  on  the  RI IS  are  orthogonal.  Taking  the  norm  of  (5.1  I )  and  using  the  fact  that  (c)  is 
perpendicular  to  Pba,  one  obtains  (cross  terms  are  0) 


4"  =  Hr+r*4. 


(5.12) 


which  is  the  Pythagorean  theorem,  for  this  right  triangle  in  function  space. 
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Figure  5.1.  Pythagorean  theorem  in  L  ?  function  space. 


As  a  brief  aside,  note  that  the  Schwarz  inequality  merely  recognizes  that 
say  that 


in  (5.12),  to 


r 


J 


df 


(5.13) 


so  that  one  has 


2>  |<  a,b  >1 


with  equality  iff  [e  =  0] — i.e.,  when  (a)  is  parallel  to  (b). 


(5.14) 
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For  application  here,  transpose  (5.12)  to  obtain 


(5.15) 


which  can  be  written  out  as 


df\ 


a(f)~intrbU 


dj  . 


(5.16) 


The  term  on  the  RHS  of  (5. 16)  is  (>()),  with  equality  iff: 

a(f)  =  yb(f);(Vf),  for  some  y.  (5.17) 

(Recall  that  sets  of  measure  0  arc  ignored  here.) 

One  may  now  apply  (5.16)  to  (5.5).  multiple  times.  Studying  the  form  of  (5.5),  one  finds  six  terms 
on  the  RHS,  each  being  a  product  of  two  O  factors.  To  write  subsets  of  terms  from  the  RIIS  of  (5.5)  in  the 
form  of  the  LHS  of  (5.16),  one  should  use  two  terms  per  subset.  One  needs  to  choose  the  correct  three 
pairs  of  terms,  and  to  choose  a(f)  and  b(f)  appropriately  for  each  pair.  Hypothesize  that  a(f)  and  b(l)  both 
consist  of  products  of  S]  (  / )  and  S2  (  /  ) ,  and  their  conjugates.  Note  from  the  LHS  of  (5.1 6)  that  in  both 
terms  there  are  second-order  products  both  in  a(f)  and  in  b(f).  The  hypothesis  then  forces  one  to  group 
terms  from  the  RHS  of  (5.5)  having  identical  total  counts  of  indices  1  and  2.  This  yields  the  following 
grouping: 


A  v’g.{det[//(^)]}  =  [OllllO,l2l <1>I12I] 

(f> 

+  [(I).2.2<I)2222  -^22.2^1222] 
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One  still  must  identify  a(f)  and  b(f)  for  each  bracketed  term.  Recall  from  (5.8)  that  (b)  must  have  positive 
norm.  One  arrives  at  the  following  identifications,  which  end  up  working  out  well,  as  will  be  seen. 

First  lenn  on  RHS  of  (5.18) 


6(/)  =  S1'(/)S,(/);<I(/)  =  S1'(/)S1(/). 


From  (3.29,  5.18,  5.19)  and  the  RHS  of  (5.16), 


A,  = 


"(/)-%#*(/) 


df 


=  [J|S,(/)|4<//] 


5;(/)5,(/)-<^rl,y>i5l(/)i 


df 


One  obtains  from  (5.17)  that  (  A,  >  0) ,  and 

a,  =  o  iff  s'2  (f)st  (/)  =  rX (f)st  (/);(v/),(F/, ) 
[s;(/)-r, s;(/)]s,(/)=o 

S2  (/)  =  y'\$\  {  f  )\  whenever  5,  (/)  *  0,(F/,). 


Second  term  on  RHS  of  (5.18) 


*(/)  =  *?  (./>(/)  =  S](f). 


From  (5. 18,  19)  and  the  RHS  of  (5. 16)  (as  above). 


(5.19) 


(5.20) 


(5.21) 


(5.22) 
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(5.23) 


A2  =[||5,(/)|V 


15,-1 


df 


One  obtains  from  (5.17)  that  (A,  >0),  and 


A:  -0  ///522(/)  =  r25,2(/);(y/),(^,). 


Third  tenn  on  RHS  of(5.18) 


/.(/)  =  S,(/)5;(/);a(/)  =  S,-(/)5!(/), 


Analogously  to  above  is 


Aj  = 


\\sAftdf 


From  (5.17),  (A3  >0),  and 


a,  =o  iff  s;(f)s2(f)=rMf)sl(fY(vf),(vr>), 
[AV(/)-rJs;(/)]s.(/)= o 

5,  (/)  =  rlsi  (f);  whenever  S2  (/)*  0,(F>3). 


From  (5.18,  21, 24,  27),  one  concludes  that 


/I vg.jdet [//(<*>)])  =A,  +A2+A,  >0, 


(5.24) 


(5.25) 


(5.26) 


(5.27) 


(5.28) 
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with  equality  iff 


■%(/)  =  rS,(f);(v/),(vr).  (5.29) 

One  simple  way  to  obtain  (5.29)  is  to  observe  that  (5.24)  shows  that  £,(/)  and  52(/)  must  have 
identical  nonzero  supports,  and  then  (5.21. 27)  easily  lead  to  (5.29). 

Now,  one  may  apply  (5.28)  to  (3.25),  to  reach  conclusions  for  the  MIMO  clutter  covariance.  First, 
recall  that 


zlvg.{det[//(^)]}=zlvg.jn  (5.30) 

=^.{a,[//(^)]a2[//(^)]},  <5-3') 


because  T  -  2  in  this  section.  But  then,  from  (3.25),  (5.28,  31), 


stvg.{A,  [l(*)].l,  [!((*)]}  =  R-  Avg.U  [//(*)],!,[//(*)]} 

<t>  <j> 

-  R2  z)yg.|det[// (^)]} 

<t> 

—  R  ( Aj  +  At  +  A3 )  ^  0, 


with  equality  iff  (5.29)  holds.  The  {  A. }  for  the  sum  on  the  RIIS  of  (5.32)  are  found  in  (5.20,  23,  26). 

Note  that  the  mathematical  expressions  in  (5.20,  26)  indicate  that  the  penalty  for  differences 
between  Sj  (  / )  and  S2  (/)  is  particularly  severe  near  spectral  peaks  of  either  waveform,  which  tends  to 
magnify  the  differences.  This  would  suggest  that  waveforms  having  flat  spectra  might  be  advantageous. 
Also,  note  that  A2  can  be  kept  low  by  just  keeping  Sj2(/)  and  S*  (/)  proportional,  which  might  be 
advantageous,  if  feasible. 
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The  result  in  (5.32)  ean  be  converted  into  lower  bounds  on  f,  itself"  or  on  (/i,  'A, ).  using  the 
expression  for  H  [<p)  in  (5.2),  and  recalling  that 


=  +  [//(«()].  <5.33) 

The  trace  is  periodic  in^ ,  but  its  average  value  over  (p  or  its  maximum  value  over  tp  are  easy  to  obtain 
from  (5.2). 

Because  £(^)  has  only  nonnegative  eigenvalues,  the  average  of  f  {<p)A^  {(/))  in  (5.32)  can  be 
0  iff  4  [<p)A^  {(p)  =  0  for  all  <j).  But  l(^)  is  of  rank  1  or  0  iff  \  (<p)  =  0.  Thus,  the  result  in 

(5.32)  implies  the  following  statement,  which  is  utilized  and  extended  in  the  next  section: 


For  T  =  2,  the  MIMO  clutter  covariance  X(^)  is  of  rank  1  or  0  for  all 
<p  iff  the  waveforms  are  proportional.  (5.34) 

The  next  section  extends  this  statement  to  T  >  2. 
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6.  IMPACT  OF  UNEQUAL  WAVEFORMS  ON  THE  SECOND  EIGENVALUE 
OF  THE  MIMO  CLUTTER  COVARIANCE,  WHEN  T  >  2 


The  purpose  of  this  section  is  to  extend  (5.34)  to  T  >2.  The  proof  involves  two  steps.  The  first 
step  establishes  the  equivalence  of  the  three  mathematical  conditions  shown  in  Figure  6.1,  regarding  the 
MIMO  clutter  covariance  and  the  waveform  cross-correlation  matrix.  The  second  step  uses  the 
equivalence  of  the  three  conditions  to  reach  the  desired  result. 

The  equivalence  of  Conditions  (A-C)  will  be  established  by  proving,  in  order,  the  four  implications 
labeled  in  Figure  6. 1 . 


MIMO  Clutter  Covariance  has 
Rank  1  or  0 


Condition  A 


Rank  E(z  ZH)]  =  1  or  0;  for  all  4> 


Condition  C 


Implication  1 

Condition  B 


Delay-Similar  Waveform  Cross¬ 
correlation  Matrix 


Delay-Similar  Vandermonde 
Mapping 


C(r)  =  C(0)r(r);  \r 

for  all  x,  for  some  scalar 
function  y(r);  [ y(0)  =  1] 


C(r)v  A</>)  =  x(<l>)ri{r  ,</>): 

for  all  (t,4>).  for  some  x(^), 
and  some  scalar  function 


0,0)  =  1] 


Figure  6. 1.  Equivalence  of  conditions. 


Proof  of  Implication  I 


Condition  B  states  that 


C(0iV(^)  =  *(^)'7(^):V(r,^), 


(6.1) 


where 
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/7(0,^)  =  1  . 


(6.2) 


Note  that  (6. 1 , 2)  imply  that 


x{t)  =  C(O)vT(0). 


(6.3) 


Also,  note  that  x(^)  may  be  0  for  some  values  of  ()). 
From  (3.21, 22)  then. 


H {4)  =  *(0)[j]/7(r’^)|  ((/)), 


(6.4) 


so  H  [(/))  has  exactly  one  nonzero  eigenvalue,  if  is  not  0  (recall  [6.2]),  and  H  {(/))  is  a  0  matrix  if 

X\(t>)  is  0.  By  (3.25)  then,  X(^)  has  less  than  two  nonzero  eigenvalues,  and  thus  X(^)  has  rank  1  or 
rank  0,  establishing  Condition  A. 

Proof  of  Implication  2 

Assume  that  Condition  B  is  not  true.  It  shall  be  shown  that  then  Condition  A  is  also  not  true.  Note 
that  G(r,^)  in  (3.22)  is  positive-semidefinite  for  all  (r,0).But  G(r,^)  is  the  kernel  of  the  integral  for 
II  in  (3.21).  If  Condition  B  is  not  true,  then  instead  of 


one  would  have,  for  some  <f>0 , 


C(r)vr(^)  =  £(r,4)7( r,4). 


(6.5) 


with  x(r,^n)  and  T)  ( T, )  guaranteed  to  be  nonzero  for  some  values  of  (t)  and  with  x (r,^»0) 
guaranteed  to  change  its  direction  as  (r)  varies.  But  then,  in  (3.21,  22), 


H  ( 4a )  =  J * (  ) \n ( r,  )|  xH  (r, </>0 ) dr 


(6.6) 
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The  integrand  is  positive-scmidefinite,  Hennitian;  hence,  the  integral  can  grow  only  in  the  ordering  of 
positive-semidefinite  matrices  over  the  integration  range  of  (r) .  Because  x(r,^n)  changes  direction  as 
r  varies,  H  [<f>  0 )  must  have  rank  greater  than  1 .  and  (3.25)  ensures  that  I  ( </)„ )  also  has  rank  greater  than 
I .  Thus,  Condition  A  is  not  met. 

Proof  of  Implication  3 

If  Condition  C  is  met.  then 


C(r)  =  C(0)7(r);(Vr), 


(6.7) 


r(o)  =  i. 


(6.8) 


Then 


C(r)v7.(^)  =  C(0)vr(^)y(r);(Vr). 


(6.9) 


To  satisfy  Condition  B,  define 


*(0)  =  C(O)vr(^),  (6.10) 

(6.ii) 


Proof  oflmplication  4 

Condition  B  gives  (6.1,  2).  It  will  be  shown  that  this  implies  Condition  C,  which  is  (6.7.  8).  First,  by 
(3.31),  the  autocorrelation  functions  icnn  (r)}  satisfy 

c_(0)>0;(»  =  1 . T).  (0  12) 

In  (0.3),  x(4>)  was  identified  as 

*(«>)  =  C(0)vr(«I).  (0  13) 
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One  must  keep  in  mind  that  x{(/))  may  be  0  ,  for  some  values  of  <(>.  Pick  an  arbitrary  value,  r0 ,  for  x.  It 
will  be  shown  that  (6.7)  holds  for  r0 .  Then,  (6.8)  is  a  simple  consequence  of  (6.7),  proving  Condition  C, 
because  r0  was  arbitrary.  Substituting  (6.13)  in  (6.1),  Condition  B  implies  that 


C(ro)vr(^)  =  C(0)vr(^)77(r0,^);(V<z>).  (6.14) 


Because  ^  is  periodic,  one  can  write  r/ (r0,^)  as  a  Fourier  series: 


m-  -oo 


(6.15) 


Substituting  (6. 1 5)  into  (6. 14),  one  obtains 


C(r„)&(*)  =  [c(0)vr(*)] 


Z  1-(T o)e 


jm<t> 


(6.16) 


Consider  the  nlh  row  of  (6.16).  On  the  LHS,  one  has  a  sum  of  complex  exponentials,  having  forms  in  the 
set 


je/0,e/*,...,e'(r (6.17) 

Because  (6.14)  is  assumed  to  be  true,  the  coefficients  j  Tjm  ( r0 ) j  must  result  only  in  frequencies  on  the 
RHS  of  (6. 16)  that  are  included  in  the  set  (6. 1 7).  This  must  hold  for  every  row  of  (6.16)  (i.c.,  for  all  n). 

Consider  the  first  factor  in  brackets  on  the  RHS  of  (6.16).  In  the  nth  row  of  this  vector,  the 
coefficient  of  exp[/(«-l)^]  is  cm)(0), which  is  positive,  per  (6.12).  Thus,  if  one  chooses  (>?  =  l), 
the  coefficient  of  exp[/0^]  is  c,,(0).  In  order  not  to  generate  frequencies  outside  of  (6.17),  it  is 
necessary  that 


(ro)  =  0  ;  w  <  0. 


(6.18) 


Now,  let  (n  =  T).  One  concludes  similarly  that 
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Hm  (ro)  =  0;w>0. 


(6.19) 


Thus,  only  //0  ( r„ )  can  be  nonzero,  so  from  (6. 1 6) 

C(r0)vr(^)  =  C(0)vr(^)70(r0);(V^).  (6.20) 

This  must  be  true  for  all  r  ,  because  r„  was  arbitrary. 

Finally,  one  may  pick  )  to  yield  a  Vandermonde  basis  for  C7  : 

[M$) . Y.t{<Pt)\-v-  <6-2l> 

Then,  by  (6.20) 

C(r)K  =  C(0)^o(r);(Vr),  ,6.221 

and  one  can  postmultiply  by  V"1  to  obtain 


C(r)  =  C(0)<fc(r);(Vr).  (6.23) 

To  obtain  (6.7),  one  merely  identifies 

r{r)  =  nM-  (6-24> 

Recall  that  (6.7)  implies  (6.8),  so  Condition  C  has  been  proven,  when  Condition  B  was  assumed.  Thus, 

implications  1-4  have  been  established.  The  equivalence  of  Conditions  A-C  follows.  In  an  aside,  the  fact 

that  Condition  A  implies  Condition  C  is  another  way  of  looking  at  why  uncorrelated  signals  cannot 
support  good  clutter  cancellation.  As  noted  in  (3.37).  uncorrelated  signals  also  have  very  mismatched 
autocorrelation  function  shapes,  making  their  waveform  correlation  matrices  not  delay  similar.  Thus, 
Condition  C  docs  not  hold,  and  hence  Condition  A  cannot  hold.  Also,  the  fact  that  Condition  C  implies 
Condition  A  tells  why  ideal  waveforms  would  have  yielded  a  rank  1  or  rank  0  MIMO  clutter  covariance, 
if  such  waveforms  had  existed  over  extended  range  domains.  This  explains  the  good  S1NR  performance 
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shown  in  Section  1,  for  simulated  ideal  waveforms.  In  fact,  when  ideal  waveforms  were  simulated,  C(0) 
was  forced  to  equal  1T,  so  that  in  Condition  C 


C(r)=/rr(r);(Vr), 


and  thus  in  (3.22), 


G(r,^)  =  vr(^)|/(r)|2  v"  (</>). 


Then  in  (3.21), 


^(^)  =  Vr(^)[{lr(r)|2  d*\ Yt  (^)» 


and  in  (3.20),  using  (3.19), 

S(^)  =  5(^)vr(^)[J|r(r)|2  rfr]v? (*)*"(*) 

=[M^)®7r]M^)[jMr)l2  ^]yt 


(6.25) 


(6.26) 


(6.27) 


(6.28) 


Thus,  the  MIMO  clutter  covariance  for  ideal  waveforms  involves  the  Vandermonde  structure  lor  the 
(single)  steering  vector  of  the  clutter  ray.  This  further  explains  the  narrowness  of  the  notch  in  the  SINR 
performance  as  the  velocity  of  the  target  changes,  making  (()  in  (6.28)  vary  away  from  ^  in  (3.15),  for  the 
(Vandermonde)  target  signature.  This  is  the  same  sharpness  that  would  be  achieved  for  a  SIMO  system 
using  real  receivers  at  all  the  locations  of  the  MIMO  virtual  array  elements,  as  discussed  in  the  text  for 
Figure  1.4.  Note  that  the  structure  of  in  (6.28)  effectively  would  be  the  same  whether  the 

transmitters  or  receivers  were  spaced  close  together.  This  is  why  the  SINR  performances  of  these  three 
designs  are  the  same,  as  mentioned  in  the  text  discussion  of  Figure  1.4. 

It  is  now  shown  that  (5.34)  can  be  extended  to  T  >  2.  First,  if  the  T  waveforms  are  all  proportional, 
then  all  the  ct/  (r)  are  all  proportional  autocorrelation  functions.  Trivially,  Condition  C  is  met,  and  hence 
Condition  A  is  met,  by  the  equivalence  of  conditions.  Next,  if  any  pair  of  the  T  waveforms  are  not 
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proportional,  then  one  knows  that  if  these  were  the  only  two  transmitted  waveforms,  the  MIMO  clutter 
covariance  would  not  be  of  rank  1  for  all  <|>,  by  (5.34)  .  But  then  Condition  A  of  Figure  6.1  would  not  hold 
for  T  =  2.  Consequently,  Condition  C  would  not  hold  for  T  =  2,  by  the  equivalence  of  conditions. 
However,  the  waveform  correlation  matrix  for  T  -  2  is  a  subset  of  that  for  the  full  T;  hence.  Condition  C 
would  fail  for  the  full  T.  Thus,  Condition  A  would  fail  for  the  full  T.  So,  one  can  extend  (5.34)  to  say  that 

For  any  T,  the  MIMO  clutter  covariance  Z(^)  is  of  rank  1  or  0 


for  all  <(>  iff  all  T  waveforms  are  proportional. 


(6.29) 


This  extends  the  result  of  (5.34)  to  T  >  2  .  However,  using  the  properties  of  analytic  functions  it  is 
possible  to  strengthen  (6.29)  (and  (5.34)  by  implication)  in  terms  of  their  applicability  over  the  domain  of 
(f) .  First,  one  writes  (6.29)  as  a  pair  of  statements 


and 


For  any  T,  if  all  T  waveforms  arc  proportional,  then  the  MIMO  clutter 
covariance  X(^)  is  of  rank  1  or  0  for  all  (f) .  (6.30a) 


For  any  T,  if  any  pair  of  the  T  waveforms  are  not  proportional. 

then  the  MIMO  clutter  covariance  Z(^)  must  have  rank  >1  for 

some  value  of  (j)  .  (6.30b) 

Below,  (6.30b)  will  be  strengthened  to  read: 


Note  that  the  potential  increase  in  spacing  above  a/2  for  two  arbitrary  transmitters  has  no  impact  on  the  relevant 
proofs,  because  the  proofs  dealt  directly  with  (|)  from  (3.1). 
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For  any  T,  if  any  pair  of  the  T  waveforms  are  not  proportional,  then  the 
MIMO  clutter  covariance  Z(^)  can  have  rank  1  or  0  only  for 
isolated  values  of  (f) . 


(6.31) 


The  following  argument  was  obtained  from  Keith  Forsythe  of  Lincoln  Laboratory  [19].  It  is  first 
noted  that  all  the  elements  of  the£(^>  )  are  analytic  in  ,  where  <f)i  is  obtained  by  adding  an  imaginary 
part  to^  ,  which  is  a  real  variable.  This  fact  follows  immediately,  because  these  elements  are  all  just  sums 
of  complex  exponentials  in  <f>r ,  as  shown  by  (3. 14,  19,  20,  42). 

Next,  it  is  noted  that  a  complex  matrix  W  can  be  Hermitian  and  of  rank  1  or  0  iff 


W  =  wa  w11  ;(for  some  complex  vector  vv  and  real  scalar  a).  (6.32) 

This  is  true  via  eigendecomposition.  It  is  shown  next  that  (6.32)  is  equivalent  to  the  two  following  sets  of 
conditions  on  W. 


=0  y(jXt,m).  (6.34) 

The  functions  {/?,  * j  m  (•)}  are  analytic  in  their  arguments.  The  (-)Jare  not  analytic  functions,  but 
that  won’t  matter  below.  The  fact  that  (6.32)  implies  (6.33,  34)  is  easy  to  show.  First,  (6.33)  just  says  that 
W  is  Hermitian,  which  is  implied  by  (6.32).  Second,  (6.32)  implies  that 


=(»<,  w>,»;  a! 

=  0;(V j,k,l,m). 


(6.35) 


proving  (6.34).  Next,  it  is  shown  that  (6.33,  34)  imply  (6.32).  The  proof  requires  a  few  steps.  Consider 
any  complex  matrix  W  satisfying  (6.33,  34).  If  W  =  0,  then  (6.32)  is  met  when  a  -  0.  So  say  that  W  *  0, 
so  a  *  0.  If  the  elements  of  the  nth  column  of  W  are  all  0,  then  (6.33)  implies  that  the  elements  of  the  n,h 
row  of  W  are  also  all  0.  The  union  of  that  row  and  column,  which  has  the  shape  of  the  crossing  of  a 
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horizontal  line  with  a  vertical  line  within  W,  satisfies  the  condition  in  (6.32),  with  w„  set  to  0.  Thus,  one 
may  remove  that  column  and  row  from  W,  and  consider  whether  (6.32)  is  met  for  the  remainder  of  the 
matrix.  Assume  that  W  now  has  had  all  such  0  columns  and  corresponding  0  rows  removed.  Thus,  every 
column  of  W  now  has  some  nonzero  entry. 

One  may  now  show  that  the  reduced  matrix  W  in  fact  has  no  zero  entries.  Suppose  that  element  W|k 
is  0.  One  can  then  obtain  a  contradiction.  It  is  known  that  the  k,h  column  also  contains  a  nonzero 
element — call  it  W(k.  Pick  an  arbitrary  column,  (m  *  k).  It  is  known  from  (6.34)  that 

WikWlm=WjmWlk.  (6.36) 


The  LHS  is  0  if  Wjk  -  0,  and  it  is  known  that  W,k  is  nonzero,  so  W/m  is  0.  Thus,  the  entire  j:h  row  is  0. 

By  (6.33),  the  entire  j'1'  column  is  also  0.  But  this  row  and  column  would  have  been  removed  in  the  initial 
reduction  of  W,  which  is  a  contradiction.  So,  the  reduced  W  has  no  zero  entries. 

Because  the  reduced  W  has  no  zero  entries,  one  may  now  consider  arbitrary  (  j,k,Lm)  and 
examine  the  implications  of  (6.36).  If  one  divides  both  sides  by  Wfk  W,m  (which  cannot  equal  0),  one 
obtains 


W  k 

_ M 

w 

n  Ik 


w 


w 


I  m 


(6.37) 


As  this  equality  ranges  over  j  and  f  ,  it  shows  that  the  kth  and  m,h  columns  are  proportional.  Because  k 
and  m  were  arbitrary,  all  columns  are  proportional,  and  W  is  of  rank  1.  But  (6.33)  said  that  W  was 
Hcrmitian,  so  it  must  have  the  form  in  (6.32)  for  its  eigendecomposition.  This  establishes  the  equivalence 
of  (6.32)  to  (6.33,34). 

One  may  now  return  to  strengthening  (6.30b).  When  any  pair  of  the  T  waveforms  are  not 
proportional.  (6.30b)  says  that  X(^)  must  have  rank  >1  for  some  value  of  <j).  Consider  such  a  value  of  c(). 
Then  1 {<(>)  cannot  be  written  as  in  (6.32),  and  hence  (6.33,  34)  cannot  all  be  true.  However,  I(^)  is 
Hcrmitian,  so  (6.33)  is  true.  Thus,  (6.34)  must  be  violated.  In  other  words. 

for  some{j,kJ,m).  (6.38) 


Extending  <)>  to  complex  </>t ,  (6.38)  says  that  for  some  tf>i  on  the  real  axis 
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h/.kj.n,  [I(^, )]  *  0;  for  some  (j,k,£,m). 


(6.39) 


But  the  composite  functions  of  </>t  in  (6.39)  are  all  analytic  in  <j>c ,  because  all  elements  of  )  are 
analytic  in  </>t  .  Thus,  if  any  of  these  functions  is  not  zero  for  a  point  in  the  complex  plane,  then  it  can  only 
be  0  at  isolated  points  in  the  plane  [26,  pp.  739-743],  hence  only  at  isolated  points  on  the  real  axis  (which 
means  real  <(>  values).  Because  of  the  equivalence  of  (6.33,  34)  to  (6.32),  one  sees  that  (6.32)  can  be  true 
only  for  for  isolated  values  of  (j).  But  (6.32)  is  equivalent  to  being  Hermitian,  of  rank  1  or  0. 

However,  it  is  already  known  that  l(^z ))  is  Hermitian  (3.20),  so  X((Z))  can  only  be  of  rank  1  or  0  for 
isolated  values  of  <|).  Condition  (6.31)  follows,  as  desired. 
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7.  PROSPECTS  FOR  IMPROVING  THE  CANCELLATION  OF  CLUTTER 

OVER  EXTENDED  DOMAINS 


This  section  briefly  discusses  three  different  potential  techniques  for  improving  the  cancellation  of 
clutter  over  extended  domains.  The  techniques  arc  adaptively  weighted  slow-time  taps  (i.e.,  Space-Time 
Adaptive  Processing,  or  STAP  [18]),  adaptively  weighted  fast-time  taps,  and  nonconventional  reception 
(that  is,  reception  not  assuming  the  use  of  matched  filters  for  the  T  waveforms).  The  discussion  focuses 
on  the  case  of  uncorrelated  waveforms,  as  considered  in  Section  4.  A  more  comprehensive  treatment  of 
this  subject  would  consider  the  joint  problem  of  waveform  design  together  with  the  processing  options 
listed  above. 

STAP  is  examined  first,  for  the  case  of  repetitive,  pulsed-Doppler  waveforms.  The  SINR 
performance  examples  given  in  the  text  utilized  two  pulse-repetition-interval  (PRI)  staggers  [18,  p.  125], 
One  could  consider  adding  more  staggers,  or  utilizing  fully  adaptive  STAP.  In  the  latter  case,  for  the 
hypothesized  range  cell  of  the  target,  the  sampled  matched-filter  outputs  lie  in  the  (.1-0)  domain  of 
(antenna  number.  MF  number,  pulse  number).  Consider  the  following  question.  For  uncorrelated  pulse 
wavelorms,  will  any  of  the  outputs  in  the  2-D  domain  of  (antenna  number,  pulse  number)  tor  MF#k  have 
any  correlation  with  any  of  the  outputs  for  MF#m,  for  k  *  m?  The  answer  is  no.  as  is  now  discussed.  In 
Section  4,  the  discussion  considered  only  clutter  on  the  target  isoDop.  Flowever.  note  that  I  !(<}>)  in  (4.4)  is 
a  scaled  identity  matrix  for  all  values  of  <j>.  Ilence,  for  each  pulse  alone,  with  no  preDoppler  filtering,  one 
can  conclude  (similarly  to  the  discussion  at  the  end  of  Section  4)  that  the  processed  antenna  outputs  for 
MF#k  are  all  uncorrelated  with  those  for  MF#m,  unless  k  =  m.  For  short  coherent  processing  intervals, 
the  cross-correlation  gain  patterns — see  (3.17) — on  the  ground  corresponding  to  the  various  combinations 
of  transmit  waveforms  and  MF  numbers  do  not  change  significantly  for  different  pulses.  Thus,  it  is  clear 
that  there  w  ill  be  little  correlation  between  any  of  the  outputs  in  the  2-D  domain  of  (antenna  number, 
pulse  number)  for  MF#k  with  any  of  those  for  MF#m,  unless  k  -  m.  Hence,  the  2-D  domain  of  outputs  for 
each  matched-filter  number  must  be  weighted  to  form  its  own  null,  which  effectively  increases  the 
required  number  of  degree  of  freedom  for  clutter  cancellation.  Consequently,  fully  adaptive  STAP  is  not 
expected  to  restore  ideal  SINR  performance. 

Adaptive  weighting  of  fast-time  taps  at  the  MF  outputs  is  now  examined.  In  this  case,  the  MF 
outputs  would  no  longer  be  sampled  only  at  the  hypothesized  target  range,  but  also  at  other  ranges  in  the 
vicinity  of  the  hypothesized  target  range,  for  adaptive  weighting.  The  adaptive  weighting  would  be  over 
the  domain  of  (antenna  number,  MF  number,  range  sample).  The  concept  of  the  argument  here  is  similar 
to  that  for  slow-time  taps.  For  uncorrelated  waveforms,  the  question  is  whether  any  of  the  outputs  in  the 
(2-D)  domain  of  (antenna  number,  range  sample)  for  MF#m  will  have  any  correlation  w  ith  any  of  the 
outputs  for  MF#n,  when  m  *  n.  The  answer  is  no,  as  shown  below.  Return  to  considering  only  the  clutter 
on  the  isoDop  to  the  hypothesized  target.  For  any  specific  range  sample,  including  the  target  range,  the 
argument  at  the  end  of  Section  4  shows  that  the  antenna  outputs  for  MF#m  are  all  uncorrelated  w  ith  those 
for  MF#n,  unless  m  =  n.  However,  one  must  now  consider  the  correlations  between  samples  at  different 
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ranges  and  MF  outputs.  This  requires  some  modifieations  to  the  previous  derivation.  The  MF  outputs  at 
the  target  range  were  expressed  in  (3.17),  which  is  repeated  here. 


i  =  (7.1) 

The  discussion  in  Section  3  considered  the  covariance  of  z.  Now  one  must  examine  the  cross-correlation 
of  z  and  zE,  where  e  is  a  lime  offset  corresponding  to  sampling  the  MFs  at  a  range  somewhat  offset  from 
that  of  the  target.  The  impact  of  the  sample-time  offset  is  to  offset  sk(  )  in  (3.6),  resulting  in  an  offset  of 
the  gain  pattern,  C(x),  being  put  on  the  clutter  on  the  isoDop.  Thus,  similar  to  (7.1), 

It  =[E«(^)®/r]JC(r  +  f)T/  (</>)u(r)dr.  (7.2) 

One  must  examine  the  cross-correlation  between  z  and  z t.  : 

I  ,(*)^  £•[«?].  <7-3> 

Analogous  to  (3.20-22),  one  obtains  here 

!,(,/>)=  (7.4) 

where 

=  (7.5) 

and 

G£(^,f)  =  C(r)vr(^)vr  (t  +  s).  (7.6) 

The  key  step  is  to  modify  the  expressions  in  (3.28,  29),  to  account  for  the  sample-time  offset: 
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Jc,y(  T)cu(T  +  e)dT 

=  I [/■ S>  (/) SJ  ( f) exP I  J2nf r) #] [{ sl ( r ) S'  (r ) exp { j 2ny ( r  +  e ) j dy ]  dr 


=  If5*  (/)5/  (/)5*  O')5*  (r )[J exp [j2x(f-y)  r}  dr 
=  J 5*  (. f)sj  (. /)5*  (f)s;  (/) exp { -jlnfe } df 


exp  j  -jlnye }  df  dy 


(7.7) 


=<VM- 


(7.81 


One  obtains,  analogously  to  (3.41), 

Zc«/(r)c-(r+^)exp{M/-*)}.  <7-9) 

*=i  i=i 


and  from  (7.5,  8)  one  has 


[H-:  =  X  i>xp{M''-*)K™,*  M- 

*=i  /=i  (7.10) 

But  for  uncorrelated  waveforms,  (3.36)  showed  that  the  waveforms  share  no  frequency  support.  Thus,  by 
(7.7)  one  has  analogously  to  (4.1)  that 


(*)  =  0,  unless m  =  i  =  n  =  k. 


(7.11) 


Thus,  similarly  to  (4.2), 


[>,((*)]  = 

m  =  n 

L  6  x' 

V  0;  otherwise  ^ 

(7.12) 


II  one  substitutes  (7.12)  into  (7.4)  and  uses  (3.10),  one  finds  that 
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M  -  [v«  ( <f>)  ®  Zr  ]  (  T <D_  ( * ) )  [ v«  ( </> )  ®  Z r  ] 

=  [v,(^)®/r][i®^g(ro_(^))][v"(<Z))®/r] 

=  [v*  (<%"  (^)] ® diag {TQ>mmmm (e )). 


This  expression  implies  that  none  of  the  antenna  output  samples  for  MF#m  correlate  with  any  for  MF#n, 
whatever  the  time  shift  between  samples.  Thus,  cancellation  of  the  clutter  must  use  weights  that  cancel 
the  antenna  outputs  from  MF#m  (with  some  benefit  from  extra  range  samples,  perhaps)  and  also  cancel 
the  outputs  from  MF#n.  This  requires  extra  degrees  of  freedom  and  indicates  that  adaptive  weighting  of 
extra  fast-time  samples  should  not  be  expected  to  restore  ideal  performance. 

Consideration  is  now  given  to  nonconventional  receiver  processing,  where  this  term  refers  to 
processing  not  assuming  the  use  of  matched  filtering.  Consider  the  case  of  a  pulsed-Doppler  waveform. 
The  clutter  data  for  each  antenna  element  can  be  likened  to  the  data  measured  for  SAR  imaging,  and  can 
be  considered  to  lie  in  the  2-D  (fast  frequency,  pulse  number)  domain,  with  its  2-D  Fourier  transform 
being  the  (range,  Doppler)  SAR  image  domain.  Consider  the  data  for  one  antenna  in  the  (fast  frequency, 
pulse  number)  domain.  For  the  case  of  uncorrelated  pulse  waveforms,  it  is  known  that  the  waveforms  do 
not  share  spectral  support  in  the  fast-frequency  domain  (see  (3.36)).  For  the  mlh  transmitter,  the  measured 
data  in  the  frequency-support  region  has  the  frequency  response  Sm(f)  (here  referring  to  a  single  pulse) 
imposed  by  the  transmitter.  Imagine  that  the  data  for  each  antenna  was  reversibly  normalized, 
conceptually,  by  multiplying  it  by  [s'm  (/)/ 1  Sm  (/)  |2  J,  and  that  this  was  done  separately  for  each 
transmitter’s  support  region  in  frequency.  If  the  data  were  Doppler  filtered  to  the  Doppler  frequency  of 
the  hypothesized  target,  each  antenna  element  would  just  have  data  in  the  fast-frequency  domain 
corresponding  to  clutter  on  the  isoDop  to  the  target.  The  overall  data  would  thus  lie  in  the  2-D  (fast 
frequency,  antenna  element)  domain.  Because  the  isoDop  is  just  a  ray  in  cone  angle  off  the  platform 
motion  direction,  it  can  clearly  be  cancelled  by  weighting  the  antenna-element  domain  samples  at  each 
frequency,  producing  nulls  at  each  frequency,  along  the  isoDop.  However,  to  get  cancellation  by 
combining  data  from  different  frequencies  is  another  matter.  Because  the  clutter  can  be  modeled  as  white 
in  reflectivity,  the  clutter  along  a  long  ray  in  range  on  the  ground  induces  (through  the  Fourier  transform) 
a  correlation  function  that  is  very  narrow  in  the  transform  domain  of  fast  frequency.  Hence,  samples  for 
different  transmitters  (at  different  frequencies)  will  not  correlate  well.  Thus,  clutter  cancellation  will  not 
be  able  to  rely  on  combining  data  from  different  transmitters.  This  suggests  that  the  nonconventional 
receiver  approach  also  fails  to  address  the  basic  problem  for  uncorrelatcd  pulse  signals,  while  it  still  may 
improve  performance  somewhat. 
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8.  BENEFITS  OF  A  LIMITED  RANGE-DOPPLER  DOMAIN  OF  CLUTTER 


As  discussed  in  Section  2,  the  results  in  Sections  3-7  were  generated  under  the  assumption  of  an 
extended  range  domain  of  clutter.  Of  course,  the  earth  is  of  finite  extent,  and  the  antenna-element 
footprints  on  the  ground  may  be  of  much  lesser  extent,  both  in  range  and  azimuth.  Further,  for  monostatic 
radars  the  return  power  for  a  constant  scattering  area  on  the  ground  falls  off  as  +  /-)  ' .  where  ( +  r) 
is  the  slant  range  of  the  scattering  region — recall  that  r  was  defined  in  (3.3)  as  the  offset  of  the  range  of  a 
clutter  patch  from  that  of  the  hypothesized  target  at  range  /;, .  Under  some  circumstances,  these  factors 
dramatically  impact  the  achievable  SINR  performance,  for  appropriately  chosen  waveforms.  It  is  well 
known  that  waveform  sets  such  as  time-division  multiple  access  (TDMA)  are  able  to  yield  near-ideal 
performance  with  clutter  having  a  suitably  limited  range  domain.  Here,  it  will  be  shown  how  this  fact 
relates  to  the  analysis  of  the  extended-domain  case  performed  in  the  previous  sections.  The  discussion 
will  treat  only  TDMA,  although  a  number  of  different  types  of  waveform  sets  are  also  able  to  exploit 
finite  domains  of  clutter  [3,  1 5]. 

The  discussion  here  assumes  a  basic  pulsed-Doppler  radar  waveform  design.  Fach  transmitter  emits 
a  pulse  of  duration  Tp  seconds,  every  Tr  seconds.  The  transmitters  are  interleaved  in  lime,  as  shown  in 
Figure  8.1.  Thus,  for  T  transmitters  there  is  a  time  spacing  of  (Tr/T)  between  pulse  start-times.  All 
pulses  sent  by  all  transmitters  during  the  coherent  processing  interval  are  identical  in  modulation.  The 
autocorrelation  functions  of  all  T  waveforms  arc  thus  identical,  also  shown  in  Figure  8.1.  However,  the 
cross-correlation  functions  are  quite  dissimilar,  as  also  shown  there.  The  key  fact  is  that  within  the  time 
delay  of  [(Tr  /  T)  -  Tp  J ,  the  cross-correlation  functions  arc  all  0;  hence  the  waveform  correlation  matrix 
C(t)  is  delay  similar,  over  the  interval 


t\<{T,IT)-Tp. 


(8.1) 


This  says  that  Condition  C  of  Section  6  holds  over  this  interval  of  x.  Going  back  to  the  derivation  of  the 
equivalence  of  Conditions  A-C  in  Section  6,  one  may  observe  that  those  derivations  applied  for  clutter 
having  any  domain  of  x.  Thus,  Condition  A  also  holds,  so  that  the  MIMO  covariance  is  of  rank  I  or  0  for 
all  values  of  (|>.  provided  that  the  clutter  returns  are  within  the  domain  in  (8.1).  If  the  clutter  extends 
beyond  the  domain  in  (8.1),  then  the  rank  grows  beyond  I.  However,  if  the  clutter  outside  the  domain  is 
weak  enough  (c.g.,  due  to  (r0  +  /-)  losses  or  a  limited  antenna  elevation  footprint),  then  the  extended 
clutter  will  not  be  of  operational  significance,  being  too  weak  to  require  good  cancellation. 
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Figure  8. 1.  TOMA  waveforms  and  correlation  functions. 


Briefly  stated,  the  system  design  tradeoff  to  use  TDMA  to  maintain  a  low-rank  MIMO  clutter 
covariance  would  be  as  follows.  The  azimuthal  coverage  region  of  the  antennas  determines  the  angular 
domain  over  which  clutter  returns  need  to  be  cancelled.  Call  this  ±0C  in  off-broadside  angle.  The 
Doppler  cycle  shift  pulse-to-pulse  for  a  given  receiver  element  (i.e.,  the  normalized  Doppler  frequency 
[18,  p.  13])  is  given  by 


uj  =  2v,  (sin  0)7],  /A,  (8.2) 

where  va  is  the  velocity  of  the  platform,  and  0  is  the  off-broadside  angle  of  the  clutter.  To  avoid  Doppler 
ambiguities  in  the  coverage  region,  |  nr  |  must  be  kept  below  Vi.  Thus, 
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(8.3) 


2vu (sin (?  )Tr/A<  1/2  , 

T  <_J__ 

4vu(sin(9)' 

(The  reason  that  Doppler  ambiguities  must  be  avoided  is  that  these  would  add  rays  of  clutter  to  cancel, 
increasing  the  rank  ol  the  MIMO  clutter  covariance.)  To  have  a  low-rank  MIMO  covariance  matrix  using 
TDMA,  the  elevation  pattern  and  (r0  +  r)  falloff  must  be  such  that  the  clutter  requiring  cancellation  is 
limited  in  delay,  as  given  in  (8.1).  with  Tr  from  (8.3).  Thus,  the  applicability  of  TDMA  for  maintaining  a 
low-rank  MIMO  clutter  covariance  depends  on  the  system  design  in  terms  of  va,/t,  the  azimuth  and 
elevation  coverage  region,  the  number  of  transmitters,  and  the  duty  cycle. 

It  is  useful  to  note  that  the  TDMA  waveforms  are  specifically  chosen  to  have  matching 
autocorrelation  functions,  as  required  to  meet  Condition  C  of  Section  6,  and  hence  to  produce  a  rank  I  or 
rank  0  MIMO  clutter  covariance,  via  equivalent  Condition  A  of  Section  6.  The  fact  that  the  clutter  is 
range-Doppler  limited  is  not  suflicient  to  obtain  good  clutter  cancellation,  unless  the  waveforms  are 
specifically  chosen  to  exploit  that.  For  example,  fully  uncorrelated  waveforms  have  orthogonal 
correlation  functions  (3.37),  which  often  differ  markedly  in  their  primary  lobe,  in  the  immediate  vicinity 
of  the  target.  This  could  cause  increased  rank  of  the  MIMO  covariance  even  for  range-Doppler  limited 
clutter.  This  argues  against  the  use  of  waveforms  sets  (like  the  Wolf  Whistles  used  in  Figure  1 .6)  that  are 
designed  for  low  cross-correlations  at  all  delays.  As  the  domain  of  the  clutter  region  increases,  the  best 
design  of  MIMO  waveforms  for  minimal  degradation  of  clutter  cancellation  while  achieving  other 
beneficial  properties  remains  an  open  issue.  Airborne  experimental  data  will  be  used  to  assess  the 
resulting  performance. 
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9.  SUMMARY 


The  use  of  uncorrelated  waveforms  for  MIMO  GMTI  is  desirable,  because  it  diffuses  the  transmit 
energy  over  a  broad  azimuthal  coverage  region,  leading  to  increased  effective  coherent  integration  lime. 
Unfortunately,  as  has  been  observed  in  computer  simulations  and  now  analyzed  herein,  when  the  clutter 
occupies  an  extended  domain  of  delay,  the  clutter  cancellation  achieved  by  using  practical,  low- 
correlation  sets  of  waveforms  is  significantly  worse  than  that  for  ideal  uncorrelated  waveforms  (i.e., 
uncorrelated  waveforms  having  matching  autocorrelation  functions).  Such  ideal  sets  of  waveforms  do  not 
exist  over  extended  domains  of  delay  (they  do  exist  for  suitably  limited  domains  of  delay).  For  extended 
domains  of  delay,  actual  uncorrclated  waveforms  must  have  nonoverlapping  frequency  support  (3.36), 
giving  them  orthogonal  autocorrelation  functions  (3.37).  The  clutter  on  the  ground  is  weighted  by  these 
autocorrelation  functions,  causing  a  loss  of  correlation  in  the  MIMO  clutter  covariance.  It  is  proven  in 
Section  4  that  for  T  uncorrelated  transmitters  the  MIMO  clutter  covariance  matrix  has  T  equal,  positive 
eigenvalues,  causing  clutter  cancellation  problems.  (For  a  single  transmitter,  the  rank  of  the  comparable 
clutter  covariance  matrix  is  1  or  0 — see  3.39.)  The  outputs  of  the  matched  filters  corresponding  to  each 
transmitter  are  uncorrelated  with  the  outputs  for  all  other  transmitters,  as  seen  at  the  various  receiver 
outputs.  Thus,  the  clutter  cancellation  is  performed  only  by  the  receivers  (Section  4).  For  a  critically 
spaced  virtual  array,  if  the  receivers  arc  closer  together  than  the  transmitters,  then  the  notch  in  the  SINR 
loss  versus  hypothesized  target  velocity  is  wider  than  the  ideal  ease.  If  the  receivers  arc  farther  apart, 
multiple  narrow  notches  are  experienced,  because  of  array  ambiguities.  In  effect,  the  dependence  of 
clutter  nulling  on  the  receive  array  design  would  thus  be  the  same  for  MIMO  as  for  SIMO. 

If  the  waveforms  are  allowed  to  be  somewhat  correlated,  it  was  shown  in  Section  6  that  if  any  pair 
of  the  T  transmitters  are  not  proportional,  then  the  rank  of  the  MIMO  clutter  covariance  must  be  greater 
than  I  for  all  but  isolated  values  of  the  direction  of  the  clutter-ray  isoDoppler  to  the  hypothesi/ed  target. 
For  the  special  case  of  T  =  2,  the  average  value  of  the  product  of  the  top  two  eigenvalues  of  the  MIMO 
clutter  covariance  was  shown  (in  Section  5)  to  have  a  quadratic  penalty  whenever  the  waveforms  are  not 
proportional.  Unfortunately,  the  penalty  emphasizes  differences  between  waveforms  at  the  frequencies 
where  either  waveform  is  most  powerful.  This  suggests  the  use  of  flat-spectrum  waveforms. 

A  brief  discussion  was  given  in  Section  7  of  the  potential  benefits  of  adaptively  weighting  either 
slow-time  taps  (i.e.,  STAR)  or  fast-time  taps.  Also,  the  potential  use  of  receivers  that  do  not  employ 
matched  fillers  for  the  waveforms  was  considered.  It  was  found  that  the  fundamental  problem  caused  by 
the  use  of  uncorrelated  signals  is  not  addressed  by  any  of  these  techniques,  although  some  performance 
improvements  might  be  obtained. 

For  limited  domains  of  clutter  in  range  Doppler,  the  example  of  TDM  A  waveforms  was  considered, 
showing  how  they  could  yield  ideal  performance  in  the  context  of  the  analysis  done  herein.  The  system 
parameter  requirements  for  ideal  performance  with  this  waveform  type  were  discussed.  The  antenna 
coverage  footprint  in  azimuth  and  range  on  the  ground  are  important,  as  well  as  the  platform  v  elocity,  the 


center  frequency,  and  the  duly  cycle  (Section  8).  As  the  domain  of  the  clutter  grows,  the  best  design  of 
MIMO  waveforms  for  minimal  degradation  of  clutter  cancellation  and  other  beneficial  properties  remains 
an  open  issue.  Airborne  experimental  data  will  be  used  to  assess  the  resulting  performance. 
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